NUMBER VARIANCE OF RANDOM ZEROS ON COMPLEX MANIFOLDS 



BERNARD SHIFFMAN AND STEVE ZELDITCH 

Abstract. We show that the variance of the number of simultaneous zeros of m i.i.d. Gaussian 
random polynomials of degree N in an open set U C C m with smooth boundary is asymptotic 
to TV'™ -1 / 2 v mm Vol(dU), where v mm is a universal constant depending only on the dimension 
m. We also give formulas for the variance of the volume of the set of simultaneous zeros in 
U of k < m random degree- A polynomials on C m . Our results hold more generally for the 
simultaneous zeros of random holomorphic sections of the 7V-th power of any positive line 
bundle over any m-dimensional compact Kahler manifold. 
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1. Introduction 

This article is concerned with the asymptotic statistics of the number Af^(p^, . . . ,p%) of 
zeros in an open set U C C m of a full system {p^} of m Gaussian random polynomials (or 
more generally, of sections of a holomorphic line bundle over a Kahler manifold M m ) as the 
degree N — > oo. In earlier work [SZlJ, we proved that the zeros become uniformly distributed 
in U with respect to the natural volume form. The main result of this article (Theorem 11.11) 
gives an asymptotic formula for the variance of A/^(pf, . . . for open sets with piecewise 
smooth boundary. We also give analogous results for the volume of the simultaneous zero set 
of k < m polynomials or sections (Theorem 11.41) . Our results show that the zeros of a random 
system are close to the expected distribution, i.e. number statistics are 'self-averaging', and 
moreover the degree of self-averaging increases with the dimension. 



Research of the first author partially supported by NSF grants DMS-0100474 and DMS-0600982; research of 
the second author partially supported by NSF grants DMS-0302518 and DMS-0603850. 
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To introduce our results, let us start with the case of polynomials in m variables. By homog- 
enizing, we may identify the space of polynomials of degree N in m variables with the space 
H (CF m ,O(N)) of holomorphic sections of the iV-th power of the hyperplane section bundle 
over CP m . This space carries a natural SU(m + l)-invariant inner product and associated 
Gaussian measure 7^. To each m-tuple of degree iV polynomials (p± , . . . ,p^), we associate its 
zero set {Pi(z) = ■ ■ ■ = p^(z) = 0}, which is almost always discrete, and thus obtain a random 
point process on CP m . We denote by 



Z p f,..., p N :={zeM: p» \z) = ■■■= p»{z) = 0} 



the set of zeros and by 



the sum of point masses at the joint zeros. It easily follows from the SU(m + l)-invariance of 
7jv that the expected value of this measure is a multiple of the Fubini-Study volume form, i.e. 



E 



N m [-u FS ) , (2) 



, 11 

where u;fs = f <9<91og \z\ 2 is the Fubini-Study metric on CP m . Here, E(X) denotes the expected 
value of a random variable X. 

Given a measurable set U, the random variable counting the number of zeros of the polyno- 
mial system in U is defined by 

<(pf , . . . ,pj[) := #{* G 17 : p?{z) = • • • = rffc) = 0}. 

Clearly, is discontinuous along the set of polynomials having a zero on the boundary dU . 
Integrating (J2J) over U gives 

E{M$) = N m Ul^Y ■ (3) 

Formula (ED has a counterpart for Gaussian random holomorphic sections of powers of any 
Hermitian holomorphic line bundle (L, h) with positive curvature G/j over any m-dimensional 
compact Kahler manifold M with Kahler form u = §©ft. The Hermitian inner product 

(a u e 2 )=f h N {a u <J 2 )—^ m , <r lt a 2 G H°(M, L N ) , (4) 



M 



induces the complex Gaussian probability measure 

d lN {s N ) = ^e-^d C) s N = J2^Sf, (5) 

on the space H°(M, L ) of holomorphic sections of L , where {S^ , . . . , S 1 ^} is an ortho normal 
basis for H°(M, L N ), and dc denotes 2d A r-dimensional Lebesgue measure. The Gaussian mea- 
sure 7at given by flUJ-flS]) is called the Hermitian Gaussian measure induced by h. The Gaussian 
ensembles (H°(M, L N ), 7^) were used in |SZ1| ISZ2| IBSZlt IBSZ2] ; for the case of polynomials 
in one variable, they are equivalent to the SU(2) ensembles studied in |BBL[ IHal INVj and 
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elsewhere. In [SZTj, we showed that the expected value of the corresponding random variable 
A/]} (where U C M) on the ensemble (H°(M, L N ) m , j™) has the asymptotics 

E«)= f (^ h y+ O(l) . (6) 



Thus, zeros of Gaussian random systems of sections of L N become uniformly distributed with 
respect to the curvature volume form ^ju: m , as iV — > 00. 

Our main result gives an asymptotic formula for the number variance 

Var(A^) :=E(A^-E(A^)) 2 

in this general setting: 

Theorem 1.1. Let (L, h) be a positive Hermitian holomorphic line bundle over a compact Tri- 
dimensional Kahler manifold M . We give H°(M, L N ) the Hermitian Gaussian measure induced 
by h and the Kahler form lu = |0/t. Let U be a domain in M with piecewise C 2 boundary and 
no cusps. Then for m independent random sections s 1 ^ G H°(M, L ) , 1 < j < m, the variance 
of the random variable 

...,&) :=#{z£U: s?(z) = ■■■ = s N m {z) = 0} 

is given by 

Var(A^) = N m ~ 1 ' 2 [u mm Yo\ 2m ^{dU)+0{N^ +£ )\ , 
where v mm is a universal positive constant. In particular, V\\ = g^//j ■ 

Here, we say that U has piecewise C k boundary without cusps if for each boundary point 
z G dll, there exists a (not necessarily convex) closed polyhedral cone K C M? m and a C k 
diffeomorphism p : V — > p(V) C M 2m , where V is a neighborhood of z , such that p(V H U) = 
p(V)(lK. By 0(iV~2+ £ ) ) we mean a term whose magnitude is less than C P N P for some constant 
C p G M + (depending on M, L, h, U as well as p), for all p > — |. 

As mentioned above, a special case of Theorem 11.41 gives statistics for the number of zeros 
of systems of polynomials of degree N. Identifying polynomials on C m of degree iV with 
H°(CF m , O(N)) endowed with the Fubini-Study metric, we obtain an orthonormal basis of 
monomials (sec [BSZ2, SZlJ): 

{G) 1/2 ^--- z t} lJm (J=(h,...j m ), \j\=h + ---+j m , (?) = (y -iA-w ) 

The polynomial case of Theorem 11.11 then takes the form: 
Corollary 1.2. Consider the Gaussian random polynomials 

p?(z 1 ,... 1 z m )= Yl CjA N j) 1/2 4 1 ■ ■ ■ (l</<m), 

{JeN™:|J|<Af} 

where the Cjj are independent complex Gaussian random variables with mean and variance 
1. Let U be a domain in C m with piecewise C 2 boundary and no cusps. Then the variance of 
the number of zeros in U of the degree-N polynomial system (pf , . . . is given by 

Var(A^) =N m ^ 2 [u mm Yo\Z\{dU) + 0{N-^) 
where Vol^,.! denotes the hypersurface volume with respect to the Fubini-Study metric on CP r ' 
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The variance Var(jV^) measures the fluctuations in the number of zeros in U of random 
systems of polynomials or sections. Theorem 11.11 implies that the number of zeros in U is 
'self- averaging' in the sense that its fluctuations are of smaller order than its typical values. 
Recalling (J6]), we have: 

Corollary 1.3. Under the hypotheses of Theorem \l.l\ or Corollary 1 1.21 

[Var(A^)l 1/2 

i yJtJj AT--4 as N -> oo . 

E(A/#) 

By Corollary O and the Borel-Cantelli Lemma applied to the sets {N- 2m (J\f^-BJ\f^) 2 > e}, 
the zeros of a random sequence {(sf , . . . , : N — 1, 2, 3, . . . } of systems almost surely become 
uniformly distributed: 

J ^tf}>( S »,...,&)-+m\Vol(U) a.s. 

Our proof of Theorem 11.11 also yields asymptotic formulas for the variance of volumes of 
simultaneous zero sets Z n n H U, where 

Z s *,..., s » :={zeM: s?(z) = ■■■ = s N k {z) = 0} . 

of k < m sections , . . . , Sf? ■ For N sufficiently large (so that L N is base point free), Z s n^ ^ s n 
is almost always a complex codimension k submanifold of M (by Bertini's Theorem) and 

Vol 2m _ 2fc (Z sf « n U) = f 1 . (7) 

1 * Jz„ N ,nu (m-k)\ 



• k 



We have the following asymptotic formula for the variance of the volume in (J7J): 

Theorem 1.4. Let 1 < k < m. With the same notation and hypotheses as in Theorem li.il 
for k independent random sections Sj 1 G H°(M, L N ) , 1 < j < k, we have 

Var(Vol 2m _ 2fc [Z sfv .. jSf f\U]) = N 2 ^ 1 ' 2 [u mk Vo\ 2m -i(dU) + 0{N^ +e ) 
where v m k is a universal positive constant; in particular, 

Vml = — g — ({m + \) . 

Theorem 11.11 is the case k = m of Theorem ll.4[ where the 0-dimensional volume is the 
counting measure: 

Voh(Z s »_ s »nU)=AfV(s?,...,s%) . 

In fact, Theorem [L4] is our means of proving Theorem ll.il as our analysis makes use of induction 
on the codimension. 

The universal constants v mk in Theorem 11.41 must be nonnegative, since variances (of non- 
constant random variables) are positive. In §4.11 we give explicit formulas for the v m k and show 
that they are strictly positive for all m and k. 

In the remainder of the introduction, we discuss related results on the variance problems 
studied in this article and indicate some key ideas in the proofs. In particular, we indicate 
which aspects of the results and methods are essentially complex analytical and which aspects 
are mainly probabilistic. 

The first results on number variance in domains appear to be due to Forrester and Honner 
|FH] for certain one-dimensional Gaussian ensembles of random complex polynomials. They 
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gave an intuitive derivation of the leading term of the asymptotic formula in dimension one, 
which is proved in Theorem 11.11 Peres- Virag have precise results on numbers of zeros of 
Gaussian random analytic functions on the unit disc for a certain ensemble with a determinantal 
zero point process |PVj . To our knowledge, there are no prior results on number or volume 
variance in higher dimensions. 

Variance asymptotics have also been studied for smooth analogues of numbers statistics, 
namely for the random variables (0Q) with ip e C°°(M). A rather simple and non-sharp estimate 
on the variance for the smooth linear statistics was given in our article [SZlJ to show that the 
co dimension-one zeros of a random sequence {s^} almost surely become uniformly distributed. 
Smooth linear statistics were then studied in depth for certain model one- dimensional Gaussian 
analytic functions by Sodin-Tsirelson |ST| ) as a key ingredient in their proof of asymptotic 
normality for linear statistics. For their model ensembles, they gave a sharp estimate for the 
variance of (Z s N,<p) and determined the leading term. (The constant j& was given for model 
ensembles in a private communication from M. Sodin.) 

We now discuss some key ideas in the proofs, and also their relation to Sodin-Tsirelson |ST] 
and to our prior work |BSZll ISZ3] . Apart from a rather routine computation (Lemma 13.31) 
for the expected value E(log |Yi| log II2I) where Yj are complex normal random variables, the 
principal ingredients in our work are purely complex analytical. The key ones are: 

• A bipotential formula for the variance current of one section (Theorem 13.11) or several 
sections (Theorem 13 .131) . and the closely related formulas for the pair correlation current 

(cf. ( !79l) ) in the codimension one case (Proposition 13. 101) ; 

• Analysis of the singularities along the diagonal of the pair correlation and variance 
currents, particularly in the point case (maximal codimension) where these currents 
contain a delta- function along the diagonal inMxM (Theorem I3.15p . This analysis 

is necessary to verify the formula K^. = [K^j (see where K^. is the pair 

correlation current for the simultaneous zeros of k random sections or polynomials of 
degree N, and to define the product [K^] A/c . 

• Application of the rapid (in fact, exponentially fast) off-diagonal asymptotics of the 
Szego kernel of [SZ2] as the degree N — > 00 to obtain asymptotics of the variance 
current and number variance. 

Let us discuss these items in more detail. The 'bipotential' for the pair correlation current 
was introduced in |BSZlj ; in the notation used here, the bipotential is a function Qn(z, w) such 
that 

A z A w Q N (z,w) = K£(z,w), (8) 
where K^ x is the 'pair correlation function' (see (1791) — ()8ip ) for the zeros of degree N polynomials 
on CP m or of holomorphic sections of L N — > M . In this article, we show that Qn{z,w) is 
actually a pluri- bipotential for the variance current Var(Z s jv) of Z s n (see Theorem 13. ip . so we 
can use Qn{z,w) to give an explicit formula (Theorems 13. lip for the variance of the zeros of 
k < m independent sections of H°(M, L N ), for any codimension k. 

The existence of the bipotential makes essential use of the complex analyticity of the poly- 
nomials and sections, in particular the Poincare-Lelong formula. It is not clear if there exists a 
useful generalization of (|S]) to the non-holomorphic setting. Moreover, the analysis of the sin- 
gularities of the variance current requires a detailed study of intersections of complex analytic 
varieties and currents of fixed degree. Part of the length of this paper is due to the lack of a 
prior reference for the relevant facts on smoothing and intersection of currents. We hope that 
the analysis in $3] will be useful in other applications. 
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The bipotential for the variance was also used implicitly by Sodin-Tsirelson [ST] , where it is 
defined as a power series in the Szego kernel for O(N) — > CP 1 . They used it to obtain the first 
sharp formula for the variance of certain model one-dimensional random analytic functions, and 
further used it to prove asymptotic normality of smooth linear statistics. 

From the bipotential formulas, and the analysis of the singularities of the variance current, 
the asymptotics of the variance are reduced to the off-diagonal asymptotics for the Szego ker- 
nel Hn(z,w) (two point function). Our main results are proved by applying the off-diagonal 
asymptotics of Hn(z,w) in |SZ2j to obtain asymptotics of Q^(z,w) and then of the variance. 
One consequence of these asymptotics, which is of independent interest, is that the normalized 
pair correlation function rapidly approaches 1 off the diagonal as N increases: 

K? m (z, w) = l + (N-°°) , for dist(*, w) > N~ 1/2+£ . (9) 

(See Corollary 13 . 1 71 for a precise statement and the definition of K^.) These asymptotics make 
essential use of the holomorphic setting, and moreover of the positive curvature of the line 
bundles involved. The two-point function Un{z,w) and consequently the correlation functions 
may behave quite differently in non-holomorphic cases (the two point kernel can decay at only 
a power law rate in real cases) or even for random polynomials of degree N on domains in 
C m with the flat metric, i.e. with an inner product independent of N. Subsequent to [SZ2J, 
sharper off-diagonal estimates for Hn(z, w) with exponentially small remainder estimates away 
from the diagonal (i.e., when dist(z, w) 3> -?=) were also given in |DLM[ [MM] , and they would 
improve the decay of correlations in (jHJ); we state the result as above, since the estimates of 
[SZ2J already suffice for our applications. 

The paper is organized as follows: In §2j we review the formulas for the expected zero 
currents and describe the asymptotics of the Szego kernel for powers of a line bundle. In £j3j we 
define the variance current (in codimension one) and introduce the pluri-bipotential Qn{z,w) 
for the variance current and study its off-diagonal asymptotics. Next, we provide our explicit 
formula for the variance (Theorem 13.111) . In §HJ this formula and the asymptotics of the pluri- 
bipotential are applied to prove Theorems 1 1 . 1 1 and 1 1 . 41 on number and volume variance. Finally, 
in the Appendix (§3]), we review and to some degree sharpen the derivation of the off-diagonal 
asymptotics in §2.21 

We end the introduction with a word on the relation of this article to its predecessors posted 
on arXiv.org. The first predecessor of this article is our preprint [SZ3j . in which we proved the 
codimension k — 1 case of Theorem 11.41 This prior article did not contain results on the point 
case in higher dimensions since, as we wrote there, "new technical ideas seem to be necessary to 
obtain limit formula for the intersections of the random zero currents Z s ." The present article 
furnishes the necessary new methods (cf. $3]). The preprint [SZ3] also extended the Sodin- 
Tsirelson asymptotic normality result for smooth statistics |ST] to general one- dimensional 
ensembles and to codimension one zero sets in higher dimensions. It remains an interesting open 
problem to generalize asymptotic normality to the point case in higher dimension. The original 
arXiv.org posting of the present article (arxiv.org/abs/math/0608743vl) also contains results 
on smooth linear statistics and on random entire functions on C m and on certain noncompact 
complete Kahler manifolds. The results from these prior preprints on asymptotic normality in 
codimension one and on smooth linear statistics are given in our article |SZ4j : the results on 
the noncompact case will be presented elsewhere. 

Acknowledgment: We thank M. Sodin for discussions of his work with B. Tsirelson. We also 
thank John Baber for computing v m 2 using formulas from the previous version of this article 
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(arxiv.org/abs/math/0608743v2), where positivity of the v m k was not established for all codi- 
mensions. Baber's computation suggested that the formulas for u m k could be simplified, and 
we were then able to revise our computations to show positivity and to obtain explicit formulas 
for the constants in all codimensions. 

2. Background 

In this section we review the basic facts about the distribution of zeros and the asymptotic 
properties of Szego kernels. 

2.1. Expected distribution of zeros. In this section, we review the formula for the expected 
simultaneous zero current of k < m independent Gaussian random sections of the tensor pow- 
ers L N = L® N of a holomorphic line bundle L over an m-dimensional complex manifold M 
(Corollary 12.31) . In order to give a simple proof of our formula by induction on the codimension 
k of the zero set, we state our result in a more general form (Proposition 12.21) . In the point 
case k = m, this formula was given by Edelman-Kostlan jEKl Theorem 8.1] (for the essentially 
equivalent case of trivial line bundles) using integral-geometric methods. 

Throughout this paper, we let (L, h) be a Hermitian holomorphic line bundle over a compact 
complex manifold M. We let S be a subspace of H°(M,L), endowed with an (arbitrary) 
Hermitian inner product. The inner product induces the complex Gaussian probability measure 

1 n 
^) = -e- |c|2 rfc, 8 = Y,CjSn ( 10 ) 

71 3=1 

on iS, where {Sj} is an orthonormal basis for S and dc is 2n-dimensional Lebesgue measure. This 
Gaussian is characterized by the property that the 2n real variables Recj, Imcj (J = 1, . . . , n) 
are independent Gaussian random variables with mean and variance h i.e., 

Ecj = 0, EcjCfc = 0, EcjCfc = 5jk ■ 

We let 

n 

U s (z,z) = E 7 (IK*)© = \\Sj(z)\\l , zeM, (11) 

3=1 

denote the 'Szego kernel' for S on the diagonal. We now consider a local holomorphic frame 
cl over a trivializing chart U, and we write Sj = jj&L over U. Any section s65 may then be 
written as 

n 

s = (c,F)e L , where F = (f u . . . , f k ) , (c, F) = ^ cjfj . (12) 

3=1 

1 

If s = feL, its Hermitian norm is given by ||s(2:)||/ l = a(z) 2 \ f(z)\ where 

a(z) = || eii (z)||fc 2 . (13) 
Recall that the curvature form of (L, h) is given locally by 

Q h = dd log a , 
and the Chern form C\(L, h) is given by 
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Using standard notation, we let S p,q (X) and V p,q (X) denote the spaces of C°° (p, g)-forms 
and compactly supported C°° (p, g)-forms, respectively, on a complex manifold X, and we let 
V'P' q (X) = V m -P' m - q (Xy denote the space of (p, g)-currents on X; (T, ip) = T(ip) denotes 
the pairing of T E V' p ' q (X) and ip E T> m ~ p ' m ~ q (X). If Y is a complex submanifold of X 
of codimension p, we let [Y] E V P,P (X) denote the current of integration over Y given by 
([Y], ip) = f Y ip. (For a further description of currents on complex manifolds, see |GHt Ch. 3].) 

We now suppose that S is base point free; i.e., the set {z E M : s(z) = 0, Vs E S} is 
empty. By Bertini's theorem (or by an application of Sard's theorem), for almost all fc-tuples 
(si, . . . , Sfe) G S k , the simultaneous zero set {z E M : Si(z) = ■ • • = Sk(z) = 0} is a complex 
submanifold of M of codimension k, and we let ^ lv .. jSfc := [{si(2;) = • • • = Sk(z) = 0}] denote 
the current of integration over the zero set: 



(Z Sl _ SkJ ip)= / <p(z), ip E X> m ~ ' m ~ (M) . 

J { Sl (z)=-=s k (z)=0} 

The current Z sx ^ Sk is well-defined for almost all Si, . . . , and we regard Z su _ >Sk as a current- 
valued random variable. For the point case k = m, Z slj _ >Sm is a measure-valued random 
variable: 

(Z Sl _ Sm ,ip)= <peV(M). 

S 1 (z) = --- = S m (z)=0 

The current of integration Z s over the zeros of one section s E S, written locally as s = fe^, 
is then given by the Poincare-Lelong formula (see |GH[ p. 388]): 

Z s = —dd log |/| = ^Eldd log + Cl (L,/0 , (15) 

7T 7T 

where the second equality is a consequence of (fl3l)-( fl4l) . 

We begin with a general form of the probabilistic Poincare-Lelong formula from [SZlj (see 
also [BSZll IBSZ2] ) for the expected value of a random zero divisor: 

PROPOSITION 2.1. Let (L,h) be a Hermitian line bundle on a compact Kahler manifold M. 
Let S be a base-point-free subspace of H°(M, L) endowed with a Hermitian inner product and 
we let 7 be the induced Gaussian probability measure on S. Then the expected zero current of 
a random section s E S is given by 

E 7 (Z S ) = ^dd\ogU s (z,z) + ci(L,h) . 

The formula of the proposition was essentially given in [SZ1[ Prop. 3.1]. For completeness, 
we give a short proof: It suffices to verify the identity over a trivializing neighborhood U. As 
above, we let {Sj} be an orthonormal basis for S, we write Sj = fjeL (over U), and we let 
F = (/i, . . . , /„). As in [SZlj . we then write F(z) = \F(z)\u(z) so that \u\ = 1 and 

log|< C ,F)| = log|F| + log|<c,«>|. (16) 



Thus by (JTS1) . we have 

(E 7 (Z s ),ip) = [ (log I (c,F) lddip)d 1 



-1 



7T 

for all test forms ip E P n *- 1 '' n ^ 1 (t7) 



(\og\F\,ddip) + J ( l °S\{c,u)\,ddip) d-y 



NUMBER VARIANCE OF RANDOM ZEROS ON COMPLEX MANIFOLDS 



A key point is that (c,u(z)) is a standard (mean 0, variance 1) complex Gaussian random 
variable for all z G U (since u(z) is unit vector) and hence Eflog |(c, u(z))\) is a universal 
constant C independent of z. 



Thus 



(log|(c, u)\,dd<p) 



dl{c) \ log \{c,u)\ddip 



M 



M 



c 



log|(c,u)|d7(c) 
dd(p = 0. 



dd ip 



Fubini's Theorem can be applied above since 



'MxS 

Therefore 



/ \log\(c,u)\dd<p\ d lN {c)= ( [ \log\(\\±e~^d() (I \ddy\) < +oo . 

JMxS \JC n J \JM J 



E 7 (Z sy 



2tt 



-<9<91og|F| 2 



2tt 



dd (log^H^II^ + AMoga) 



(17) 
□ 



The formula of the proposition then follows from ([T5]) and pT|) . 

Remark: Proposition 12.11 also holds if M is noncompact and without the assumption that S 
is base point free or even finite dimensional. 

Next, we give a general result on the expected value of the simultaneous zero current of k 
independent random holomorphic sections: 

Proposition 2.2. Let (L, h) — > M, S C H°(M, L), and 7 be given as in Proposition ^. 1\ and 
let 1 < k < m. Then the expected value of the simultaneous zero current of k independent 
random sections S\, . . . ,Sk in S is given by 

E 7fe (Z Sl ,..., Sfe ) = h^ddlogU s (z,z) + c l (L,h) S j . 

The proposition is a formal consequence of Proposition 12.11 and the independence of the 
sections Sj, but needs a proof since the wedge products of currents is not always defined. We 
give here a simple induction proof without using the theory of wedge products of singular 
currents. 



,m—k 



(18) 



Proof. Let uo be the Kahler form on M. We first note that 

(Z Sl _ Sk ,u m - k ) = [ Cl {L,h) k Au 

J M 

whenever the Z s . are smooth and intersect transversely. The identity ffTSl) is a consequence 
of the fact that the current Z Sl) ... Sfc and the smooth form ci(L,h) k are in the same de Rham 
cohomology class. Equation (fTSj) can also be verified by induction: the case k = 1 follows 
immediately from the Poincare-Lelong formula (1151) and the fact that u is closed; assuming the 
result for k 

(z Sl n z. 



1 sections on Z Sl , we have 



S'2. 



UJ 



m—k\ 



c 1 (L,h) k - 1 Au r 



d(L,h) Ac 1 (L,h) k - 1 Au r 
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which gives ( 1181) for k sections. 

Now let if G T> m ~ k,m ~ k (M) be a test form. By (Tl8l) and the formula for the volume of complex 
submanifolds ([7]), we then have 



(Z ai ,...,s h ,<p)\ 



Z B -, ,...,s h 



< sup || || Vo\(Z Su 



for almost all si, . . . , s&. Thus the random variable (Z Slt ,,, tSk , </?) is L°°, so its expected value is 
well defined. 

We must show that 

= / « fc Ayp, (19) 



where 



a = -r-^dd log Us (z, z) + c x (L, /i) . 

Z7T 



We verify ({19]) by induction on k: For = 1, Proposition 12 . II yields ({19]) . Let k > 2 and suppose 
that (Tl9l) has been verified for k — 1 sections. Choose si E <S such that Z Sl is a submanifold, 
and let M' = Z Sl , s'j = Sj\m', <S' = S\m'- We give S 1 the Gaussian measure 7' := p*7, where 
p : S — > iS' is the restriction map, and we note that 

= Ey (ll^)liy = E 7 (\\s(z)\\l) = U s (z',z') , for z' e M' . 
By the inductive assumption applied to M',S', and noting that Z Sli ... jSfc = Z s i ...y , we have 

/ (Z ai ,...,s k ,<p) d-y(s 2 ) ■ ■ -dj(s k ) = Ey h -i(Z s > 2 s,cp) = \ 
Js*- 1 ' Jz Sl 

We average ( {20]) over si and apply Proposition 12.11 to conclude that 



c/" 1 A <p . (20) 



(^5 1 ,... )8Jk ,y)d7(«0"- d T(*fc)= / (Z si ,a k 1 A^)d7(si) = (a,a k 1 A = / or A <p, 

which gives ( 1191) . □ 

2.1.1. Powers of a positive line bundle. We now specialize Proposition ^. 2l to our case of interest. 
We suppose that the line bundle (L, h) is positive, i.e. the (1, l)-form c\(L,h) is everywhere 
positive definite, and we give M the Kahler form uj = = irci(L,h). Recall that it is 
a consequence of the Kodaira embedding theorem that for sufficiently large integers N, the 
spaces of global sections of the tensor powers L N = L® N of the line bundle are base point free. 
(In fact, the global sections give a projective embedding [GH] §1.4]. The Kodaira embedding 
theorem is also a consequence of the Tian-Yau-Zelditch theorem [Cal ITI] IZe2j ; see ({291.) We 
recall that the Hermitian metric h on L induces Hermitian metrics h on L , and we have 

Cl (L N ,h N ) = N Cl (L,h) = —w. 

We give H°(M, L N ) the Hermitian inner product induced by the metrics h, u, as defined by 
(TJJ; this inner product induces the Hermitian Gaussian measure 77V given by (jSJ). Considering 
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the spaces Sn = H°(M, L N ), we have the Szego kernels (on the diagonal) 

d N 

U N (z,z) := U Sn (z,z) = \\Sf(z)\\ 2 hN , (21) 

i=i 

where {S^,...,S^} is an orthonormal basis for H°(M,L N ) with respect to the Hermitian 
inner product (@J. These Szego kernels were analyzed in |Ze2l IBSZl} \SZ1\ ISZ2] by viewing 
them as orthogonal projectors on C 2 (X), where X — > M is the circle bundle of unit vectors of 
L~ x . We give this description of fLy in §2.21 below. 

Applying Proposition 12.21 to the line bundles L N and the spaces H°(M,L N ) of holomorphic 
sections, we obtain: 

Corollary 2.3. Let (L,h) — » (M,cj) be as in Theorem \ and let 7jv be the Hermitian 
Gaussian measure on H°(M, L ). Then for 1 < k < m and N sufficiently large, we have 

^ddlogU N (z,z) + —uj . 

2.2. Off-diagonal asymptotics for the Szego kernel. As in |Ze2[ISZl] IBSZl] and elsewhere, 
we analyze the Szego kernel for H°(M, L N ) by lifting it to the circle bundle X— >M of unit vectors 
in the dual bundle L^ 1 — > M with respect to h. In the standard way (loc. cit.), sections of L N 
lift to equivariant functions on X. Then s 6 H°(M, L N ) lifts to a CR holomorphic functions 
on X satisfying s(e z9 x) = e tN0 s(x). We denote the space of such functions by 7i 2 N {X). The 
Szego projector is the orthogonal projector 11^ : C 2 {X) — ► TC%(X), which is given by the Szego 
kernel 



U N (x,y) =Y,Sj r {x)S?{y) (x,y e X) . 

i=i 

(Here, the functions are the lifts to TC%(X) of the orthonormal sections S^; they provide 
an orthonormal basis for 7i 2 N (X).) 

Further, the covariant derivative Vs of a section s lifts to the horizontal derivative V^s of 
its equivariant lift s to X; the horizontal derivative is of the form 



7)1 



V>-^g-A§>, (22) 

For further discussion and details on lifting sections, we refer to [SZlJ. 

Our pluri-bipotential for the variance described in §31 is based on the normalized Szego kernels 

P N ( Z , W ):= |n ", (z '"' )l , , (23) 

U N (z,z)2U N (w,w) 



2 



where we write 

\U N (z,w)\ := |IIjv(a;,y)| , z = ir(x), w = n(y) G M . 

In particular, on the diagonal we have Hn(z, z) = Hn(x, x) > 0. Note that EEjv^, z) = U.s(z, z) 
as defined in (EI} with S = H°{M, L N ). 

In this section, we use the off-diagonal asymptotics for H N (x,y) from |SZ2] to provide the 
off-diagonal estimates for the normalized Szego kernel Pn(z,w) that we need for our variance 
formulas. Our estimates are of two types: (1) 'near-diagonal' asymptotics (Propositions 12.71 - 
12. 8p for Pn(z,w) where the distance dist(z, w) between z and w satisfies an upper bound 
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1 /2 

dist(z,w) < b (— jr-} (b G (2) 'far-off-diagonal' asymptotics (Proposition 12.61) where 

dist(^^)>6(M) 1/2 . 

To describe the scaling asymptotics for the Szego kernel at a point zq G M, we choose a 
neighborhood U of zq, a local normal coordinate chart p : U,z — > C m , centered at zo> and a 
preferred local frame at z 0; which we defined in |SZ2j to be a local frame such that 



IM«)||/i = l-il|pWII^ + -- • (24) 
For ii = (ui, . . . , u m ) G p(U), 9 G (— 7r, 7r), we let 

p(ui, . . • , M m , 0) = e^P" 1 ^)) e x > (25) 

|e L (p 1 (u))\h 

so that (ui, • • • , u m , 9) G C m x R give local coordinates on X. As in |SZ2j . we write 

U z ^(u,9;v,ip) = IL N {p{u,6),p{v,cp)) . 

Note that depends on the choice of coordinates and frame; we shall assume that we are given 
normal coordinates and local frames for each point zq G M and that these normal coordinates 
and local frames are smooth functions of z . The scaling asymptotics of n^(w, 9; v, ip) lead to 
the model Heisenberg Szego kernel 

n*(z, 9- w, <p) = e m ^ y S k (z)S^) = ^ime-rt+N^-i ( W »+m«) (26 ) 

fceN m 

of level N for the Bargmann-Fock space of functions on C m (see [BSZ2J). 

We shall apply the following (near and far) off-diagonal asymptotics from [SZ2J: 

Theorem 2.4. Let (L,h) — > (M,u) be as in Theorem [7771 and let z G M. Then using the 
above notation, 

'A AI- m Y\ z ° ( ® • v XL} 

= Uf(u, 9; v, if) [l + YLx N- r / 2 p r (u, v) + N-^^R^u, v) 
where the p r are polynomials in (u, v) of degree < 5r (of the same parity as r ), and 

\V j R Nk (u,v)\ < C jk£b N £ for \u\ + \v\ < by/ log N , 
for e, b G M + , j, k > 0. Furthermore, the constant Cjkeb can be chosen independently of 

Zq. 



ii) For b > \fj + 2k + 2m , j,k > 0, we have 

\VITIn(z,w) \ = 0(N~ h ) uniformly for dist(z, iu) > b „ ^, 

Here V j R = { Qu k>q v k" '■ \ K '\ + \ K "\ = i}> and V{ = (V h ) j denotes the j-th iterated hori- 
zontal covariant derivative; see (j22j) . Theorem 12.41 is equivalent to equations (95)-(96) in [SZ2J, 
where the result was shown to hold for almost-complex symplectic manifolds. (The remainder 
in (i) was given for v — 0, but the proof holds without any change for v ^ 0. Also the state- 
ment of the result was divided into the two cases where the scaled distance is less or more, 
respectively, than N 1 ^ 6 instead of \AogN in the above formulation, which is more useful for 
our purposes.) A description of the polynomials p r in part (i) is given in [SZ2] , but we only 
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need the k = case in this paper. For the benefit of the reader, we give a proof of Theorem 
1231 in m 

Remark: The Szego kernel actually satisfies the sharper 'Agmon decay estimate' away from 
the diagonal: 

Vm N {z, 0; W,(f)=0 (e-AiVN&stiz,^ ; ■ > q ( 27 ) 

In particular, 

\U N (z,w)\ = O (e-AVNSBt( z ,w)^ _ (2g) 

A short proof of (I28I) is given in |Bel Th. 2.5]; similar estimates were established by M. Christ 
[Oh] . H. Delin [Dij, and N. Lindholm[Li]. (See also jDLMl [MM] for off-diagonal exponential 
estimates in a more general setting.) We do not need Agmon estimates for this paper; instead 
Theorem 12.41 suffices. 

It follows from Theorem I2.4( i)) with k = 1 that on the diagonal, the Szego kernel is of the 
form 

U N (z,z) = ^N m (l + 0(N~ 1 )), (29) 

which comprises the leading terms of the Tian-Yau-Zelditch asymptotic expansion of the Szego 
kernel |Cal [Til IZe2j . Applying (1291) . we obtain the asymptotic formula from |SZlj for the 
expected simultaneous zero currents: 

Proposition 2.5. |SZll Prop. 4.4] Let (L, h) — > (M,u) be as in Theorem\TJ\ and let 1 < k < 
m. Then for independent random sections s± , . . . , in H°(M, L ), we have 

N k 

V(Z s ?,..., s «) = ^ + 0(N^). 
Proof. By ([29]), dd\ogU N (z,z) = 0(A^ _1 ). The asymptotics for E(Z 

s N ,...,s N ) then follow from 

the formula of Corollary 12.31 □ 

We now state our far-off-diagonal decay estimate for Pn(z,w), which follows immediately 
from Theorem I2.4( ii) and (T29]) . 



Proposition 2.6. Let (L,h) — > (M, u) be as in Theorem lTl\ and let P N (z,w) be the normal- 
ized Szego kernel for H°(M, L ) given by (I23p . For b > y/j + 2k, j, k > 0, we have 



/logiV 

V ] Pn{z,w) = 0{N ) uniformly for dist(z, w) > b 



N 



The normalized Szego kernel P^ also satisfies Gaussian decay estimates valid very close to 
the diagonal. To give the estimate, we write by abuse of notation, 

P„( 2o + u, Zo + „) := P K (p- W ,p-W) = ^-J^^gL-^ . (30) 

As an immediate consequence of Theorem I2.4( i). we have: 

Proposition 2.7. Let P N (z,w) be as in Proposition ^. 6] and let z e M. Forb,e > 0, j > 0, 
there is a constant Cj = Cj(M,e,b), independent of the point zq, such that 

P n(zo + ^,Zo + ^) = e-^ 2 [l + R N (u,v)] 

\V j R N (u,v)\ < CjN- 1/2+£ for \u\ + \v\ < by/log N . 
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As a corollary we have: 



Proposition 2.8. The remainder R^ in Proposition \277\ satisfies 



\R n (u,v)\<y\ u - v \ 2n ~ 1/2+£ i \VR n (u)\<C 2 \u-v\N~ 1/2+£ , for \u\ + \v\ < by / hgN. 

Proof. Since P^ (z + u, z + v) < 1 = P^ (z + u, z + u), we conclude that Rn(u,u) = 0, 
dRN\{u,u) = 0, and thus by Proposition 12. 7[ 

\VR N (u,v)\< sup \V 2 R N {u, (l-t)u + tv)\ \u-v\ < C 2 \u - v\ N~ 1/2+£ . 



0<t<l 



Similarly, 



C 

\R N (u,v)\ < \ sup \V 2 R N (u, (l-t)u + tv)\ \u-v\ 2 < -^\u-v\ 2 N~ 1/2+£ 



2 

0<t<l 



2 



□ 



3. A PLURI-BIPOTENTIAL for the variance 

Our proof of Theorems 11.11 is based on a pluri-bipotential given implicitly in |SZ1] for the 
variance current for random zeros in codimension one. More generally, for random codimension 
k zeros, we define the variance current of Z.n s n to be the current 

Var(Z, s ») :=E(Z sfi ... iS ^Z sfi ... >sf )-E(Z sfr .. iSf )BE(Z s ,... iSf ) e V' 2k > 2k (M x M). (31) 

Here we write 

S m T = irlS A n* 2 T G V" p+q (M x M) , for S G V' P {M), T G V' q {M) , 

where 7rx,7T2 : M x M — > M are the projections to the first and second factors, respectively. 
The variance for the 'smooth zero statistics' is given by: 

Var(Z sfi ... )Sf ,^) = (Var(Z sf _ sf ), <p B y?) . (32) 

Conversely, (132]) can be taken as an equivalent definition of the variance current in terms of 

Theorem 13. II below gives a pluri-bipotential for the variance current in codimension one, i.e. 
a function Q N G L X {M x M) such that 

Var(Z s iv) = (idd) z (idd) w Q N (z,w) . (33) 

To describe our pluri-bipotential Qn(z,w), we define the function 

n=l 



Alternatively, 



G( e - A ) = -— / log(l - e- 2s ) ds , A > . (35) 



A 



2tt 

The function G is a modification of the function G defined in |BSZ1] : see f )57|) . 
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Theorem 3.1. Let (L,h) -> (M,lu) be as in Theorems{lJl Let Q N : M x M -> [0, +00) 6e 
t/ie junction given by 

1_ /• p ^(^) 2 logfl - s) 

where Pn{z,w) is the normalized Szego kernel given by (1231) . T/ien 

Var(Z s iv) = (idd) z (idd) w Q N (z,w) . 

Theorem 13.11 says that 



Q N (z, w) = G(P N {z, w)) = -— / ^ rfs , (36) 



Var(Z s iv,v?) = ( - d z d z d w d w Q N , ipMtp) = / Q N (z,w) (iddcp(z)) (idd<p(w)) , (37) 

J M J M 

for test forms tp G P m ~ 1,m_1 (M). We note that Qtv is C°° off the diagonal for iV sufficiently 
large, but is only C l and not C 2 at all points on the diagonal in M x M, as the computations 
in §3. II show. Additionally, its derivatives of order < 4 are in L m ~ £ (M x M) (see Lemma [3. 7\i . 
To begin the proof of the theorem, we write 

9 N = (Sf , . . • , S^) G H°(M, L N ) dN , (38) 

where {S^} is an orthonormal basis of H°(M, L ). As in the proof of Proposition 12. 11 we write 

M*) = I ( 39 ) 

where |^jv| := ll^ll/U) 1 ^ so \ Un \ = 1- ^ or c = ( Cl > • • • ' c d N )i we write 

/ 1 \ 1 djv 



\(c,u N {z))\ = \\(c,u N {z))\\ h N . 
Lemma 3.2. 

Var(Z s iv) = — \d z d z d w d w [ log | (c, uj\r(z))| log | (c, Wat(w)) | ^(c) . 
Proof. We write sections s N G H°(M, L ) as 

s7V = E = <c, *iv> , c = (c 1; . . . , . (40) 

Writing = Fef N , where is a local nonvanishing section of L, and recalling that 

w = -Oh = -iddlog \\e L \\ h , 

we have by (fi~5l) . 

= -991og|< C)J F)| = -991og|( C ,* Ar }|--991og|| e f 7V ||, 

7T 7T 71 

= -ddlog\(c,V N )\ + —uj. (41) 

7T 7T 

Consider the random current 

Z N := -99 log I (c, *iv>| = Z s n - -u . (42) 

7T 7T 
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It follows immediately from the definition ( 13~TT) of variance currents that 
Var(Zjv) = E(Z N B Z N ) - E(Z N ) M E(Z N ) = Var(Z s 

By ( ED , we have 



E(Z N ) = -ddlog\V N \ , (43) 

7T 



whereas by f|42|) . we have 

E(^KI^) = ~ f d z d z d w d w log | (c, log | (c, ^ N (w))\ d lN (c) 

11 JC d N 



~d z 5 z d w B w I \og\(c, V N (z))\ log|(c, ty N (w))\d lN (c) . (44) 



ft JC d N 

Recalling f[3"$j) . we have 

IogK*jv(z),c)| log | c)| = log log |^iv(^) | + log |*iv(^) | \og\(c,u N (w))\ 

+ log|^Ar(w)| hg \(c,U N (z))\ 

+ log | (c, u N (w)) \ log I (c, ujv(z)) I , (45) 
which decomposes (1440 into four terms. By (j4"3~|) . the first term contributes 



log Analog \* N {w)\ = E(Z N )mE(Z N ) 



n 2 



The c-integral in the second term is independent of w and hence the second term vanishes when 
applying d w d w . The third term likewise vanishes when applying d z d z . Therefore, the fourth 
term gives the variance current Var(Z s jv). □ 



To complete the proof of Theorem I3.1[ we use the following probability lemma, which gives 
the c-integral of Lemma 13.21 

Lemma 3.3. Let (Yi, Yg) be joint complex Gaussian random variables with mean and E(|Yi| 2 ) = 
E(|Y 2 | 2 ) = 1. Then 

E(log|Yi| log|Y 2 |) =G(\E(Y l Y 2 )\) , 

where 

7 2 1 f t2 logfl - s) 
G(t) := -j - - J — ds , < t < 1 (7 = Euler's constant) . 

Proof. By replacing Y\ with e ia Yi, we can assume without loss of generality that E(YiY 2 ) > 0. 
We can write 

Y 1 = S x , 

Y 2 = (cos 0) Si + (sin#)H 2 , 

where Hi,S 2 are independent joint complex Gaussian random variables with mean and vari- 
ance 1, and cos 6* = E(YiY 2 ). Then 

EQog|Yi|log|y 2 |) = G(cos0) , (46) 

where 

G(cosfl) = -L / loglSd loglHiCos^ + Sasin^l e - (|Sl|2+|S2|2) cGi cG 2 . (47) 

Jc 2 
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The computation of G(cos9) was essentially given in [BSZ14 §4.1]. We repeat this com- 
putation here for the readers' convenience: Write 5^ = r 1 e ia , H 2 = r 2 e l<a+l ^, so that (j47j) 
becomes 

2 f°° f°° f 2n _, 2 2 . 
G{cos9) = — III r\r 2 e ^ ri+r2 ' log r 1 log \r 1 cos 9 + r2e lip sm9\d(pdr 1 dr2 ■ 
* Jo Jo Jo 

Evaluating the inner integral by Jensen's formula, we obtain 

2tt f 2tt log(rx cos 9) for r 2 sin 9 < r% cos 9 

log \r\ cos 9 + r 2 sin 9e tlf \ dip 

2n log(r 2 sin 9) for r 2 sin 9 > r x cos 9 



Hence 



/»oo /"OO 

G(cos9) = 4 / / r 1 r 2 e~^ 1+r2 - ) log r\ log max^ cos 9, r 2 sin 9) dr\ dr 2 
Jo Jo 



We make the change of variables r± = pcosip, r 2 = psin<p to get 

G(cos9) = 4 / / p 3 e~ p log(p cos ip) logmax(p cos ip cos 9, p sin (p sin 9) cos ip simp dip dp . 
Jo Jo 

Since 

logmax(pcos<y9cos6 l ,psin(^sin6') = log(p cosp cos 9) + log + (tan <p tan 9) , 
we can write G = G\ + G 2 , where 

Gi(cos^) = 4 / / (?e~ p log(p cos <p) log(p cosp cos 9) cos p sin pdpdp , (48) 
Jo </o 

G 2 (cos6') = 4 / / p 3 e~ p log(pcos<p) log(tan<ptan#) cos <p sin <p d<p dp . (49) 

Jo Jn/2~e 



From (gHD, Gi(cos0) = Ci + L7 2 logcos6>. Substituting 

cos 9 = e~ A , 

we obtain 

G x {e- X ) = C 1 - C 2 X . (50) 
We now evaluate G 2 . Since the integrand in (j49p vanishes when p = n /2 — 9, we have 



Since 



~ 2 (e~ A ) = 4 I — log tan 6 1 J / / p 3 e~ p log(p cos p>) cos p sin pdpdp . 

\d\ J J 



^G 2 (e~ x 



log tan 6> = - — log(e 2A - 1) 



rfA to 2 rfA toV 7 1 - e- 2A 

we have 

d G 2 (e- X ) = (h + h) , 



dX ^ ' 1-e 
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where 

nn/2 
p 3 e' p2 (log p) cos y? sin (p dip dp = C 3 sin 2 6» = C 3 (l - e~ 2A ) 
.72-0 

fOO /"7r/2 



Thus 



p 3 e~ p (log cos <p) cos </? sin cp dtp dp 

/2-e 

j rw/2 i /"sinS 

- / (log cos <£>) cos cpsia(fd(p — — t log t (it 

2 Jtt/2-9 2 J 

-(sin 2 fllogsin 2 #-sin 2 #) = -(1 - e~ 2A ) [log(l - e~ 2A ) - l] . 
8 8 

^G 2 (e- A ) = ± log(l - e- 2A ) + 4 C 3 - \ . (51) 



Combining (l50 l) -( l5T|) . we have 

G(e~ A ) = C 4 + C 5 A + \ [ log(l - e" 2s ) ds . (52) 

2 Jo 

By flSD, 



G(0) = [E(log 



2 / (logr)e r rdr 



1 2 2 



Substituting A = oo in (J52|) . we conclude that C 5 = and 

G(e- A ) = -L-~ log(l - e~ 2 *) rf S , (53) 

*/ A 



or equivalent ly. 

4 4 



G(t) = tAf l0g{1 ~ S) ds (0 <*<!). (54) 

□ 



Proof of Theorem \ 3. 11 Fix points 2, w G M, and let x,y & X with 7r(a;) = 2, 7r(y) = w. We 
apply Lemma I3T31 with Yi = (c,uj^(x)), Y 2 = (c,U]s[{y)). Since \(c,U]y)\ = \(c,u N )\ on, we have 

log|Fi| = log\(c,u N (z))\, log|F 2 | = log \(c,u N {w))\ . 

To determine E(Y{Y 2 ), we note that for a random s N = ^2cjS^ G H%(X), 



E ($(:r)S(y)) = £ E ( c ^) (y) = J2s?(x)S»(y) = U N (x,y) . (55) 

Since 

(C, Wjv(x)) 



j,k=l j=l 

(c,* N (x)) s N (x) 



we have by fl55|) . 
and recalling ([23 



E(yiY 2 ) 



|$jv(a;)| ^(ar, z)V2 
njv(x,?/) 



IL N (x,xy/m N (y,y)W 
E(Y 1 Y 2 )\=P N (z,w). (56) 
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Therefore, by Lemma I3~3l and ( 1561 ). 

/ log\(u N (z),c)\ log|(« J vH, C )|d 7iV ( C ) = E(log|F 1 | \og\Y 2 \)=G(P N (z,w)). 

J<£ d N 

By (EU) and (JST 



Git) = i 



7T 2 



7 2 ' 



(57) 



and hence, recalling that G o Pjy = Q N 
1 



7T 2 Jc d 



log\{u N (z),c)\ log|(itjv(tu),c)| dry N {c) = Q N (z,w) + C . (58) 
Theorem 13.11 follows by combining Lemma 13.21 and (1581) . □ 



3.1. Asymptotics of the pluri-bipotential. We now use the Szego kernel off-diagonal asymp- 
totics to describe the iV-asymptotics for the variance current Var(Z s iv) (Lemma 13. 9j) . We also 
need to know the behavior of the variance current near the diagonal. We showed in [BSZ2] (107)] 
that the co dimension-one scaling limit pair correlation grows like 

\z — w\~ 2 near the diagonal (for dimension m > 2). Our computation of the variance cur- 
rent asymptotics also gives this growth rate for the variance current (Lemma 13.71) as well as for 
its scaling limit (Lemma 13.91) . 

We begin by noting that the pluri-bipotential decays rapidly away from the diagonal: 

Lemma 3.4. For b > y/j + q + 1, j > 0, we have 



\^Q N (z,w)\ =Ol±-) , for dist(z,w) > bx/WX 



N q J v > / - 

Proof. We recall from (136]) that Qn = G o P Nl where G is analytic at (with radius of conver- 
gence 1) and G{t) = 0(t 2 ). The estimate then follows from Proposition 12.61 with k = \_^-\ ■ D 

Applying Lemma 13.41 to the pluri-bipotential formula for the variance of Theorem 13.1} we 
conclude that the variance current decays rapidly away from the diagonal. 
We next show the near-diagonal estimate: 

Lemma 3.5. For b E R + , we have 

Q N (z ,z + ^=) =G(e- |H2 ) + 0(iV-V2+ £) ) for H < 6v /f^iv. 



f N< 

Proof. Since Pn{zq, zq) = 1 and G'(t) — > oo as t — > 1, we need a short argument: let 

A N = -\ogP N . (59) 

Recalling fl3"5"|) . we write, 



1 



CO 



F(X) := C7(e- A ) = -— / log(l - e~ ls ) ds (A > 0) , (60) 

J A 

so that 

Q N = FoA N . (61) 



By Proposition 12. 8[ 



A N ( z ,z Q + — ) = \\v\ 2 + R N (y) , (62) 
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where 

R N = -\og(l + R N ) = 0(\v\ 2 N- 1/2+£ ) for \v\ < b^logN . (63) 

By QSQD, 

< -F'(X) = - JL log(l - e~») < JL (l + l og + i) . (64) 
Since ||v| 2 + R N (v) = \v\ 2 (§ + o(iV)), it follows from that 
Q JV (z ,^o + ^) = F(i|w| 2 + J R 7V (w 

= ^(|H 2 )+o( 



1 + log + rT 

\v\ 



R N (v] 



G(e-^ 2 ) + 0(N~ 1/2+£ ) , for \v\ < b^/hg~N 



We shall use the following notation: for a current TonMx M, we write 
9T = diT + 9 2 T, dx =^dzj A , 9 2 = ^dw, ^ 



□ 



9^- ' ~ z Z-^~ ■? 

where zi, . . . , z m are local coordinates on the first factor, and W\, . . . , w m are local coordinates 
on the second factor of M x M. We similarly write 

dT = B X T + B 2 T . 

In particular, we shall write 9i9i9 2 9 2 (5jv in place of d z B z d w B w QN(z, w) to avoid confusion when 
we change variables. 

Next we compute the leading term of the iV-asymptotics of B x B 2 Qn and d\Bid 2 B 2 Q ^ . We 
choose normal coordinates at a point z G M, and we recall that in terms of these coordinates, 
we have 

Qn = F o A n , A N (z,w) = y|w - z\ 2 + R N (VN z,Vn w) , (65) 
where R^ is given by fl63|) . We now write A N (z, w) ~ B N (z, w) when 



logiV 



iV 



■Ajyfou;) - B N (z,w) = O (n ^ +£ \B N (z,w)\j for \z\ + \w\ < b 
By (1651) . we have: 

B 2 Qn(z,w) w f F'(Aiv(^w))9[(«;-z)-«;], (66) 
9 i a 2 Q JV (0,«;) w -fF"(A^,«)))a[(«i-^^-]A9[()i;-2).4 (67) 
B 1 d 2 B 2 Q N (z,w) « -^F( 3 )(Ajv(z,u;))a[(w-2) -5] A9[(iD-5) ■ io] A9[(w-z) ■ w] 

-^F"(Aiv(^w)) {99[5-w] A9[(ic-z) -tw] + 9[(w-z) -5] A99M 2 } , (68) 

and hence 



d 1 d l d 2 d 2 Q N (0,w) « ^F^(Ajv(0, w)) 9(z • u;) A 9(5 ■ u>) A 9|-u;| 2 A d\w\ 
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+i£ir(3)( Ajv ( 0) «;)) [99>| 2 A 9|uf A 9|w| 2 + 9(5 ■ w) A d\w\ 2 A 99> • u>) 
+d(z ■ w) A 99(5 ■ w) A 9|u;| 2 + d(z ■ w) A 9(5 ■ w) A 99|w| 2 ] 



+^ F"(A N (0, w)) [99(5 • w) A 99(2 • u>) + 99>| 2 A 99>| 2 ] . (69) 
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Differentiating ( 1641) . we have 

F"W = ± T^TT > ( 7 °) 

and hence 

F&(\) = 0(e- 2X ) (A > 1), (71) 

for j > 0. Furthermore, by 

F'(X) = -L fog A + 17(A), r? G C°°([0,+oo)) , 



and therefore 

F^ +1 \\) = (-l) J+l (j '~7 ! A~ j + 0(1) (A > 0), 



(72) 



2tt 2 
for j > 1. 

We now use the above computation to describe the singularity of the variance current near 
the diagonal. We first recall an elementary fact: 

Lemma 3.6. Letu G C 1 ((R p \{0}) xl"), and let 1 <j <p+q. Suppose that du/dxj G L 1 (W +q ) 
and u(x) = o (\tti(x) |~ p+1 ) ; where 7Ti : MP +q — > M. p is the projection. Then the distribution 
derivative du/dxj G V(R p+q ) is given by the pointwise derivative, i.e. 

= |% V^P(R^). 

(RP\{0})xR9 OXj J(RP\{0})xR'J OXj 

Proof. Let U £ = {x G MP +q : 1 7Ti (a?) | > e}. The identity follows by integrating u by parts 
over {7 e , and noting that the boundary term 

w</? dxi • • • dxj-idxj+i ■ ■ ■ dx p+q 

goes to zero as e — > 0. □ 

Lemma 3.7. There exist a constant C m G IR + (depending only on the dimension m) and an 
integer N = N (M) G Z + such that for N > N , we have: 

i) The coefficients of the current d^d^Q n are locally bounded functions (given by pointwise 
differentiation of Qn), and we have the pointwise estimate 



\did2Q N (z,w)\<C m N for 0< \w - z\ <byj-^-. 

ii) If m > 2, the coefficients of the current dxdxdydzQN are locally L" 1 ^ 1 functions, and we 
have the estimate 



^d^Q^w^K-^- for (Xlw-zlKbJ 1 ^. 

\w — z\ A V N 

Proof. We take z = z = 0. By Propositions I2.7H2.81 we can choose N such that 



Atv(0,w) > for H<6y-^p, iV>iV . (73) 

By applying the chain rule as in . we conclude that for each N > iV , 

VQ N (z,w) = 0(\w - z\ log \w - z\), V 2 Q N (z,w) = 0(log \w - z\), 

V j Q N (z,w) = 0(\w-z\~ j+2 ) for j>3. (74) 
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Hence, the partial derivatives of Qn of order < 3 are in L\ oc) and the same holds for the fourth 
order derivatives if m > 2. By repeatedly applying Lemma [3T61 with iri(x) = w — z, we conclude 
that the currents B\B 2 Qn an d diB\d 2 B 2 QN have locally L 1 coefficients. The upper bound in 
(i) follows from (1671) . (170]) and (1731) . and hence the coefficients of B\B 2 Qn are actually in L^ c . 
The upper bound in (ii) similarly follows from (1691) . (1721) and (1731) . and hence the coefficients 
of d 1 d 1 d 2 B 2 QN are in L™" 1 . □ 

These computations show that Qn is C 1 and has vanishing first derivatives on the diagonal 
in M x M, but Qn is not C 2 along the diagonal. Lemma fBTTT i) says that B\B 2 Qn is bounded; 
however, a similar computation shows that d 2 B 2 QN{z, w) > clog |z — for a positive constant 
c. When m > 1, d\Bid 2 B 2 QN ~ |io — z|~ 2 ; but when m = 1, BiB\d 2 B 2 QN is a measure with a 
singular component along the diagonal, and off the diagonal there is cancellation in (1691) and 
diB!d 2 B 2 QN ~ |w - ^| 2 (see [BSZTT Th. 4.2]). 

Making the change of variables 

= 



in ([6"7j) and again applying (|65p and (|72[) . we obtain the following asymptotic formulas: 
Lemma 3.8. For N sufficiently large, 

d x 8 2 Q N {z Q , z + ^) = -^F"{\\v\ 2 ) B(z ■ v) A B\v\ 2 + O(N^) , 
for < \v\ < by/\ogN. 

Proof. By (pll-pSI) and ([72]) with j = 2, we have 

F" (A N (0, ^)) = F"{\\v\ 2 ) + ([§ + 0{N- 1 / 2 ^)] \v\ 2 ) ■ O {\v\ 2 N- 1 / 2 ^) 

= F"{\\v\ 2 )+0{\v\- 2 N^/ 2+£ ). (75) 

The formula follows from (1671) and (1751) . □ 
LEMMA 3.9. For iV sufficiently large, 

-d 1 B 1 d 2 B 2 Q N {z ,z + ^ w ) = NVcir z £(v) + O(\v\~ 2 N 1 / 2+£ ) for < \v\ < b^hgN , (76) 

w/jere Var^ G T ( J 2 >(M x C m ) zs ^iuen 6y 

Var^(w) := — ^ F^(||t>| 2 ) 9(2 • v) A ■ v) A <9|t>| 2 A <9|t>| 2 

_1 F^{\\v\ 2 ) [dB\z\ 2 A <9M 2 A B\v\ 2 + <9(* • v) A <9M 2 A <9<9(;z • v) 

+d(z ■ v) A Bd(z ■ v) A B\v\ 2 + <9(z ■ v) A B(z ■ v) A <9<9|-y| 2 ] 
~\ F"{\\v\ 2 ) [Bd(z ■ v) A dB(z ■ v) + dB\z\ 2 A <9<9M 2 ] . (77) 



Furthermore, 



Var ,o M _J 0(M~ 2 ) /or M>0 

Var -W-i 0(|t;| 4 e-H 2 ) /or H>1 ' ^ 



Proof. Formula (1771) follows by the same argument as in the proof of Lemma 13.81 applying (1691) 
in place of flf>7J). The estimate fl7J) follows by applying tTTT]) — fTTSj) to (J77J. □ 

The current Var^, has L\ oc coefficients if m > 2, but contains the singular term ^uj(zo)®5 Zo {y ) 
if m = 1 (see Theorem 13.151 or jBSZll Theorem 4.1]). 
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3.2. The pair correlation current. The pair correlation current gives the correlation for the 
zero densities at two points of M. It is defined to be 

:= E(Z S » H Z a * o ) E V' 2k ' 2k (M x M) . (79) 

Thus by Corollary 12.31 

Var(Z sf _ sr ) =Ka-B(Z^ > .„^)KIE(Z^ > .„ i ^) =K^- [E(Z^) HE(Z^)] fc . (80) 
As a consequence of Theorem 13.11 we have the following formula for the case k — 1 : 
Proposition 3.10. T7ie pazr correlation current in codimension one is given by 

Kg = -d 1 did 2 d 2 Q N + EZ s n B EZ s n . 
For m > 2, the coefficients of the current are in L^ 1 . 

Proof. The formula for is an immediate consequence of Theorem 13.11 and ( IHDl) . By Lemma 
I3.7( ii). the coefficients of the current didiC^c^Qjv ar e in , and hence the same holds for 
the coefficients of K^. □ 

We note that the pair correlation function K^Zjiv), which gives the probability density of 
zeros occurring at both z and w (sec [BSZ1, BSZ2]), can be obtained from the pair correlation 
current: 

(ml Lljm ^ ml = K 2 k A ^ H ^ . (gl) 

The advantage of the pair correlation current is that, because of the independence of the sf, the 
co dimension- k correlation current is the fc-th exterior power of the corresponding codimension- 
one current, i.e. 

Kg=[Kg] Ak . (82) 

This formula is analogous to the corresponding identity E(Z s n^ s n) = E(Z s jv) fc for the expected 
value in Corollary 12. 3t both formulas hold since the sj^ are independent random sections. 
However, in the case of correlation currents, the right side of (1821) is not well defined along the 
diagonal in M x M, since is singular on the diagonal. We shall show that for k < m, the 
current K^., has L 1 coefficients and (|82p holds, with [K^] Ak given by pointwise multiplication. 
However, for the point case k = m of Theorem 11.11 the pair correlation current contains 
a singular measure supported on the diagonal (see Theorem 13.151) , and the right side of (1821) 
must be interpreted as a limit of smooth currents. The singularities of necessitate a more 
complicated proof, using a smoothing method. (One can also define in an analogous way the 

n-point correlation currents K^,, which satisfy the identity K^ fe = [K^J .) 

3.3. Explicit formula for the variance. We shall write 

= T^y urn ~ k for l <k<m-l, 

$m = 1 , 

so that 

Vol 2m - 2 fc (z s n ^ s n nU^j = [z s n ) S n , xu $fc) , for 1 < k < m . 
By Corollary 12.31 the expected volume of the co dimension- /c zero current is given by 

EjVoWa^,.^ n U)) = J {EZ sN ) k A <t> k , 



13) 



*4) 



24 BERNARD SHIFFMAN AND STEVE ZELDITCH 

where 

i - N 
EZ s n = —dd log 11^, z) + — u . (85) 

In this section, we prove the following integral formula for the volume and number variance: 

Theorem 3.11. The variance in Theorem T7% is given by: 

V&r(Voh m -2k[Z s N ^N n U]) 

= ( k ) f B i § &» A (didihfaQN)*' 1 a (ez sN m EZ sN ) k ~ j a ($ fc m $ fc ), 

j=1 \3 J JdUxdU 

for N sufficiently large, where Qn is given by (1361) . EZ s n is given by (!85|) and t/ie integrands 
are in L 1 (dU x dU) . 



In particular, for the one- dimensional case k — m = 1, we have 

Var(A^)=- / B z B w Q N (z,w) . 

JdUxdU 



Theorem 13.111 follows formally from Theorem 13.11 and equations (1801) and (|82|) . To verify the 
formula rigorously, we must show that the currents B^Qn A \d1B1d2B2Q n\ 3 (1 < j < 
which are smooth forms off the diagonal in M x M, are well defined and have L 1 coefficients; 
i.e., they impart no mass to the diagonal. To do this, we use the asymptotics of Q^(z,w) as 
\z — w\ —>■ given in Lemma 13 .71 In §U we shall use the N asymptotics of Lemmas I3.8H3.9I 
together with Theorem 13.111 to prove Theorem 11.41 

Definition 3.12. We say that a current u G V p,q (X) on a Kahler manifold X is an L 1 current 
on X if its local coefficients are L 1 functions and j x \u\ dVolx < +00. (The second condition 
is redundant if X is compact.) If {«„} is a sequence of L 1 currents on X, we say that u n — > u 
in L 1 if J x \u n — u\ dVolx — * 0. 

We shall prove Theorem 13.111 by approximating xu by smooth cut-off functions, and then 
applying the following explicit formula for the variance current: 

Theorem 3.13. Let 1 < k < m. Then for N sufficiently large, 



Var(Z s iv ____ s n) = <9i<9 2 



J^(-l)^ 1 (^jd l B 2 Q N A (d l B 1 d 2 B 2 Q N ) :i 1 A (EZ s n m EZ sN ) k 3 

where the current inside the brackets is an L l current on M x M given by pointwise multipli- 
cation, Qn is given by (|3"6T) . and EZ s n is given by flHBI) . Furthermore, Var(Z s JV ) ..,, S N ) ^ s an L 1 
current on M x M if k < m — 1 . 

We need the following smoothing result for our proof of Theorem 13.131 

Proposition 3.14. Let (L, h) — >• M be as in Theorem VLl\ with dimM = m > 2. Then there 
is a positive integer N such that the following holds for all N > N : 

• for 1 < j < in — 1, K^- is an L 1 current on M x M, and is given by the pointwise 
formula Kg = (K£)'\ 
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• for all points Pq G M x M , there exist a neighborhood Q C M x M of Pq and smooth 
forms 

S e K& E £ 2 > 2 {tt) (0 < e < 1) 

such that 

i) d l (S s K%) = d 2 (S g K%) = 0; 

ii) for 1 < j < m - \, {S £ Kgy -> Kg| n *n L 1 , as £ -> 0; 

iii) for 2 < k < m, Kg A (^Kg)* -1 -> K^| n wea% ; as e -> 0. 

We postpone the proof of Proposition 13.141 to the next section, and we now prove Theorem 
13.131 assuming Proposition 13.141 The case m = 1 of Theorem 13.131 is Theorem 13. 11 So we let 
m > 2. It suffices to consider a test form ip e x> 2m ~ 2k ' 2m ~ 2k (Vl), where Q C M x M is as in 
Proposition 13. 141 By Propositions 13.10 1 and 13.141 (recalling that BxB 2 Qn has L°° coefficients by 
Lemma [3 .7p . we have 



(K£ )¥ = lim (K^ A [iSgK^]* -1 , ip) 

= IhnJ ^B^Qn ^[SeKgf- 1 Adxd 2 ip + [EZ s n M EZ s n] A ^K^]* -1 A <p} 
Si^Qtv A (-d 1 B 1 d 2 d 2 Q N + EZ s n H EZ^ f- 1 A <9i«9 2 ^ 



+ / (EZ s iv B EZ s n) A (-di<9i<9 2 <9 2 Q;v + EZ s iv m EZ sN ) k ~ l A 9? . (86) 
Jn 

Expanding the integrand and recalling (I80p . we then have 
(Var(Z s N t s n), (pj 

= (Kg, (p)-J (EZ sN B EZ sN ) k A ^ 

= J2 Q I J) / ^Qjv A (-d 1 B 1 d 2 B 2 Q N ) j - 1 A (EZ s iv IE EZ s iv) fc - J A ft&p 
+ zZ( ■ ) B^Qn Ai-drB^B^Ny- 1 A{EZ s N^EZ s N) k -i Ad^ 

= (j) / ^Qiv A {-d l B 1 d 2 B 2 Q N ) j - 1 A (EZ sN B EZ sN ) k - j A d x d 2 y . 

This is the formula of Theorem 13.131 to complete the proof of the theorem, it remains only to 
prove Proposition 13.141 

3.3.1. Off-diagonal decay of the correlation current. We now use I3.13ll3.14l to give a more ex- 
plicit formula for the pair correlation current in higher codimension and describe its off-diagonal 
asymptotics. (The results of this section are presented here for their general interest and are 
not needed for the proof of the variance formulas.) 

We note that the correlation currents K^. are smooth forms away from the diagonal in 
M x M. We now show that K^, has no mass on the diagonal for k < m, while contains a 
'delta-function' along the diagonal: 
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Theorem 3.15. Let (L, h) — > M be as in Theorem \l.l\ Then for N sufficiently large, we have: 

i) for 1 < k < m — 1, the correlation current for Kg is an L l current on M x M given by 

Kf fe = [-d 1 d 1 d 2 B 2 Q N + EZ s n B EZ s m] k , 
where Qn is given by fl36|) . and EZ s n is given by fl85|) ,- 

ii) Kg m = {-d x d x d 2 d 2 Q N + EZ sN H EZ sJ v) m | MxM ^ A + diag,(EZ.*) m , 

w/iere A = {(z, z) : z E M} is the diagonal in M x M, and diag : M —> M x M is the 
diagonal map diag(z) = (z,z). 

Proof. Part (i) is an immediate consequence of Proposition 13.101 and the first conclusion of 
Proposition 13. 141 To verify (ii), we regard the current K% m G V' 4m (M x M) (which is of order 
by its definition) as a measure on M x M. Since Qn is C°° in M x M \ A, it follows from 
Theorem ED3 (or from (JHSJ) that = {-dxd x d 2 d 2 Q N + EZ s iv B EZ s iv) m on M x M \ A. 
Hence it suffices to show that (ii) holds on A, i.e., for any Borel set A C M, 

(KL , Xdiag(A) ) = ((EZ sN ) m , Xa) ■ (87) 

To verify (J8"7j) . we note that 

( Z < ,...,< K Z s f ,...,-5; . Xdiag(A)) = # {z E A : Si(z) = ■ ■ ■ = s*{z) = 0} = (^f , xa) • 
Taking expectations and recalling Corollary 12.31 we then have 

( K 2m , Xdiag(A)) = E (z S N ^ s n M Z s N^ s n , Xdiag(A)) 

= (eZ s n_ s „, X a) = {(EZ sN ) m , xa) , 
and therefore (ii) holds on A and hence on all of M x M. □ 

Theorem 3.16. Let (L,h) — > M be as in Theorem{TJ\ and let 1 < k < m. For b > 

\fq + 2k + 3 > 3 ; the pair correlation currents satisfy the off-diagonal asymptotics 



Kg = (EZ s n M EZ sN ) k + O (N- q ) , for dist(z, w) > ^ l ° g A 



/AT 

In particular, 

Kg = (EZ s n M EZ sN f + O (N~°°) , for dist(z, iy) > N~ 1/2+e . 
Proof. By Lemma 13.41 for b > y/q + 2k + 3 , we have 

d&d&Q N = O (j^^) , for dist(z, «;) > . (88) 



59) 



Since the statement only pertains to the off-diagonal, by Theorem 13.151 and ([581 . 



K^ fc = [-d 1 d 1 d 2 d 2 QN + EZ s n M EZ gN ] k 



EZ s n m EZ s n + O 



1 



]\fq+2k-2 



for dist(z, m) > bV W 7 . By PropositionES EZ s ivKEZ^ = 0{N 2 ), and the desired asymptotics 



'N _ 

follow immediately from fl89|) . □ 
In particular, the point case of Theorem 13.161 yields the asymptotics 
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Corollary 3.17. Let (L,h) — > M be as in Theorem and define the normalized pair 
correlation function by 

K N 

jy-N _ 

2m (EZ sN M EZ sN ) m ' 
Then for b > y/q + 5, q > 1, we have 



K? m (z, w) = l + (IT*) , for distfo w) > K l ° H A 



'N 

In particular, 

K? m {z, w) = l + (N-°°) , for dist(z, w) > N" 1 ^ . 

Corollary 13 . 1 71 can also be obtained from Theorem 12.41 (or [SZ2| (95)-(96)]), using the argu- 
ment in Section 4.1 of [BSZ2J. 

3.4. Smoothing K^: Proof of Proposition 13.141 We shall use the following fact about 
averaging currents of integration over a smooth family {Y t } of submanifolds: 

Lemma 3.18. Let X and Q be complex manifolds of dimension m and n respectively, and let 
X' be a complex submanifold of X. Let y be a complex submanifold of M x Q such that the 
projections tti : y — > M and 7r 2 : y — > Cl are submersions, and let Y t = m (tt^ 1 -^}) c X , for 
t G Q. Then for all a G V n ' n (Q), we have 

• / [Y t ] a(t) G £ P,P (X), where p = codim^ = codimx^ ; 
Jten 

• for almost all t G fi, X' C\Y t is a complex submanifold of X of codimension p' := 
p + codimX', and 

[X' n Y t ) a (t) = [X'] A I [Y t j ] a(t) G V'^'(X) . 
ten Jten 

Proof. Let y = 1 (X / ) = y fl (X' x Q), and consider the commutative diagram 

y ^ y 

</ *xS \vr 2 , (90) 

x' A x n 

where 7i[ = ir\y, and i, i are the inclusions. Since 7T2 is a submersion, Y t is a smooth submanifold 
of X for all t G fi. For a test form <p G £> m -P> m -P(X), we have 

[Y t ] a(t) , <p\ = a(t) / <p = / vr^a A vr*v? = (n U 7T*a, ip) , 
ten / Jn JY t Jy 

and thus 

/ [Y t ] a(t) = n u n* 2 a , (91) 
Jten 

which is a smooth form, since the push-forward of a smooth current via a submersion is smooth. 

Since 7Ti is a submersion, y is a submanifold, and hence by Sard's theorem, the set of critical 
values of the map n' 2 '■— i^2\y has measure zero. Therefore tc'^ 1 ^) = (X' fl Y t ) x {t} is a 
submanifold for almost all t G fi. 
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Now suppose that (p G T> m ~ p '' m ~ p ' (X). Let S denote the set of critical points of the map 
n' 2 : y — > Q, and let E = ir' 2 (S fl Supp ir'*<p), which is a closed subset of measure zero in Q. Let 
p n G C°°(fl) such that p n >0,p n /* Xn^E, and Supp p n fl E = 0. As before, we have 

p n (t)a(t) / Lp = / VT2*(p„a) A vr'/V = / TX 2 (p n a) Air'*(p ■ 
n Jx'nYt Jy^s Jy 

We claim that 

r Pn (t)a(t) I <p ^ [ a(t) [ ip, (92) 
Jx'nYt Jn Jx'nY t 

Tc' 2 *(p n a) A ir'*<p / 7r£aA7rJ>, (93) 

as n — > oo, and hence 

/ a(t) [ <p= f <aA^. (94) 

To verify (|92|) - (!94|) . we first consider the case a > and <p = (3 m ~ p \ where (3 is a (compactly 
supported) semi-positive (l,l)-form on X, and apply monotone convergence to obtain (192]) - 
(|93|) for this case. Since the right side of (|94|) is given by the integral of a compactly supported 
smooth form over y', both sides of (1941) are finite, and then (|9~2"1) - (19~3I) hold in the general case 
by dominated convergence. 
Thus, 

[X' n Y t ) a(t) ,<p)= f a(t) [ cp = [ ir'*a A i$<p = ([y'\ A vr 2 *a, vr^), 
ten / Jn Jx'nYt Jy 

and therefore 

f [X' n y t ] ar(t) = 7Tx* ([y] A vr 2 *a) . (95) 
Jten 

Furthermore, since 7Ti is a submersion, 7r* is well-defined on currents and 

7Ti*([y] A 7r 2 *a) = tt 1# (7r*[X'] A 7r*a) = [X'] A 7i u n* 2 a . (96) 
The identity of the lemma follows from fl95l) - fl96|) and fl9Tl) . □ 

We now proceed to the proof of Proposition 13. 141 By abuse of notation, we let Z s n denote 
the zero set of a section s N as well as the current of integration over the zero set. We choose 
N such that if N > N Q , the zero sets Z s n, . . . , Z s n are almost always smooth and intersect 
transversely. (This holds if the Kodaira map for H°(M, L N ) is an embedding.) 

We begin by smoothing currents (locally) on M. Let a G M be arbitrary and consider a 

coordinate chart r a : V a B r := {z G C m : ||z|| < r}, with a G K C M, r(a) = 0. We let 
?7 a = T~ l (B r / 2 ). To simplify our argument below, we choose the biholomorphism r as follows: 
Embed M C CP 9 , and choose projective coordinates (Co : ' ' ' : C?) i n CP 9 such that 

• a= (1 : : ••• : 0), 

• {x G M : CjOr) = for < j < m} = 0, 

• the projection 7T : M — > CP m , x i— > (Co(^) : • • • : Cm(^)) has nonsingular Jacobian at a. 

We choose a neighborhood of a such that 7r a is injective on V' a and 7T a (V^) C C m = CP m \ 
{Co = 0}. We then choose r > such that B r C vr a (V^), and we let V a = r~ 1 (B r ) and r a = 7r a |y a . 
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The advantage of this construction is that degree bounds in M push forward under r a to 
degree bounds in CP m . In particular, if X is an algebraic hypersurface in M, then 

r a (X n V a ) C 7r a (X), deg cpm ir a (X) = deg cpq X . (97) 

(The well-known formula (1971) is easily verified by recalling that the degree of a subvariety X in 
projective space is the number of points in the intersection of X with a generic linear subspace 
of complementary dimension.) 

Let ip £ (z) = e~ 2m "ip(z/e) be an approximate identity on C m , with ip G C°°(C m ) and Supp?/> C 
B r / 2 - We consider the local smoothing operator S £ : T> n,1 (M) — > £ l,1 (U a ) given by convolution 
in the r coordinates: 

S a £ u = u * := T a [r a Mv a ) * A] G ^(C/a) , for m G D'^M) . (98) 

(Note that T a *{u\v a ) G P /1 ' 1 (5 r ), and hence its convolution with ip £ is well-defined on B r / 2 for 
0<e<l.) 

Now suppose that P = (a, b) G M x M, and let r a : \4 S r , r& : H _B r be as above, and 
let Q = U a x £/" b . We consider the approximate identity ip £ (z,w) = ip £ {z)ip £ {w) on C 2m and we 
similarly define S £ : & 2 > 2 {M) £ 2 - 2 (ft) by 

(S , £ M)=M*^:=r*[r,(M| yaXy J*^] G £ 2 ' 2 (ft) , for w g £> /2,2 (M x M) , (99) 

where r = r a x r& : V a x V& ^> B r x B r . 
Lemma 3.19. For iV > jV , 

= E(S e (Z s ,v IE Z s *)) = E(S^Z sN B G £ 2 ' 2 (ft) . 

Proof. We have 

S^g = [E(Z S * B Z sJ v)] * ^ e = E[(Z s ,v B Z^) * Vy = E(S £ (Z S * E3 Z s m)). 
Furthermore, 

E[(Z sN B Z sJ v) * ^] = E[(Z S * * ^) B (Z s ,v * V,)] = E(S £ a Z s * B 5 £ 6 ^) . 

□ 

Lemma 3.20. Let 2 < k < m, N > N . For almost all (sf , . . . , s% ) G #°(M, L N ) k , we have 

\Z 'n Kl Z„n N B ^Jvl A • ' ' A S^f^-iV ^„Jvl — > .Z.iV „JV B Z Q N „JV , 

L Sj S 1 J tL S 2 s 2 J tL s fc s fc J s l >—> s k S l '■■■' S k 

weakly in V' 2k > 2k (Q), as e —> 0. 

Proof. It suffices to consider the case where the Z s n are smooth and intersect transversely. We 

let Yj = (Z s n x Z s n) n (V^ x Vb), and we identify V a xVb with 5 r x I?, r via the biholomorphism 

t. Under this identification, Q = B r/2 x B r/2 C C 2m . For t G C 2m , let T t : C 2m -> C 2m denote 
the translation T t (w) = w + t, so that 

= / [T_ t K,-] v(t) G £ 2 ' 2 (ft) , 

where v is the Euclidean volume form on C 2m . 
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Suppose that X' is a complex submanifold of Q. We first show by induction that X' nT^l^n 
• R T tk Yk is a complex submanifold of fl for almost all t 2 , . . . ,t k , and 



\x'} a s £ y 2 a ■ • • a s e y 



[X' n T t2 Y 2 n • • • n T tfc y] ^ e (t 2 ) • • • ^ e (t fc ) i/(t 2 ) a • • • a v(t h ). (100) 

To verify ffTOQj) for fc = 2, we let y = {(z, t) G ft x : z - t G y}. Then y is smooth and the 
two projections n 1 : y — ► f2, 7r 2 : 3^ — ► ^ are submersions. Furthermore, 7r 1 (7r 2 _1 {t}) = T t Y 2 . 
Hence by Lemma 13.181 with X = Q and a = i/) e v, the intersection X' fl T t2 Y 2 is a complex 
submanifold for almost all t 2 G Q, and (11001) holds for k = 2. For the inductive step, let k > 2 
and suppose that fllOOH has been verified for k — 1. Let ii, . . . , tfc-i be parameters in fl such 
that X' fl T t2 Y 2 fl ■ ■ ■ PI T tk _ 1 Y/ C _ 1 is a complex submanifold of fl By Lemma 13.181 with X' 
replaced by X' fl T t2 Y 2 fl • • ■ fl T tkl Y k -\ and y — {(z, t) E fl xfl : z — t E Y k }, we conclude that 
X' fl T t2 Y 2 n • • • fl T tfe y, is a complex submanifold of f2 for almost all t k , and 



x' n F t ; n ■ ■ • n y*~\ n y* ^(4) = x' n y; n - • • n y*- ] 



AS £ Y k 



Integrating over ti, . . . , tfc_i and applying the inductive assumption, we obtain fllOOp . 
Setting X' = y in (PT0U1) . we have 

[y] a s £ y 2 a • • • a 5 £ y fe 



[y n T t2 y 2 n ■ • • n T th Y k ] ^ £ {t 2 ) ■ ■ ■ ^ e {t k ) u(t 2 ) a • ■ ■ a u(t k ). (101) 



Now choose Eo > such that y , T f2 y, . . . , T ffc Yfc intersect transversely whenever < e for 2 < 
j < k. Let 9? G z>*n-2fc,2m-2*(n) be a test form g ince the submanifolds y fl T t2 y 2 n • • • fl T tfe y 
vary smoothly as the parameters t 2 , . . . ,t k vary in the £o-ball, it follows from an argument using 
the implicit function theorem that the map 



(t2 



,tfc) 



([ynr^yn.-.n^y],^) 

'yinTi 2 y 2 n-nri fc Y fc 

is continuous (and in fact is C°°) for \tj\ < Eq. Therefore by fllOip . 

([y] A5 e y 2 a--- Asy^) -> ([y nyn---ny], <p) as ^ 

i.e., 

[y] a 5 £ y 2 a • • • a s^y. -»• [y n y 2 n • • • n y] weakly, as e -> o 



o 



(102) 

□ 



Lemma 3.21. There exists a positive constant C < +oo suc/i i/ia£ /or a// X > Xo, if G 
p2m-2fe,2m-2fc(^ and < e < 1, we have 



< CN 2k \\tp 



[Z s n M Z s n] A S £ [Z s n M Z s n] a 
/or a/mosi a/Z (sf , . . . s k ) G tf°(M, L JV ) fe . 

Proof. Fix si, . . . , Sfc so that Z s jv, . . . , intersect transversely, and let K- = (Z s jv x -Z^) fl 
(14 x Vj,), as in the proof of Lemma 13. 20[ By (jlOip it suffices to show that 

([y n r t2 Y 2 n • • • n T 4fc y] , < cx^imu , (103) 



NUMBER VARIANCE OF RANDOM ZEROS ON COMPLEX MANIFOLDS 31 
for almost all (t 2 , ...,£&) G Q^ 1 . To verify (11031) . it in turn suffices to show that 

VoWxm (Yx n T t2 r 2 n ■ ■ ■ n r ffc r fe ) < cw 2fc . (104) 

We write tj = (t'ptf), Z<f = Z s n n K, Z) = Z s n n H, so that T^.Y} = T^Z? x T t ,,Z). Then 

VoImxm Ci n T t2 Y 2 n • • • n r tfc r fc ) 

= VoW(z? n T t , z 2 a n • • • n 7V k ^) Voi M (z? n T q z\ n • • • n r^) . (105) 

Since r a : <^-> 7r a (Z s Ar) C CP m is injective and the translations T t '. extend to automorphisms 
of CP m , we have 

Voi M (£? nr^n-n T t < k Z a k ) < c x Voi C p™ (^(z^) n T^ a {z s n) n • • • n T^ a {z s n )) 



n deg cp m7r «( z . 



/^l A A —our" "»\- sv 

3=1 



However, by ( 1971) . 

deg cpm 7t a (Z s N) = deg C¥q Z s n = —— T [ u™~ q l = [ N a(L, h) A u™^ 1 = C 2 N, 

3 3 ^ J z N n Jm 



and hence 

Voi M [zi nr^n-n T t ,z£) < c 3 N k . (ioe) 

The bound (ITU!) follows from ffT05|) -ffT06l). □ 

Lemma 3.22. Let fj G L^ c (IR fc ) for 1 < j < n. Let /J = /j * where ipe is a compactly 
supported smooth approximate identity. Then 

n n 

lift ^ fl/i as e^O. 

j=0 j=0 

Proof. We use the generalized Holder inequality: 

n 

X)- = l =► ll/i'--/»lli<ll/ilU---|l/»lk- (107) 

3=1 Pj 

We can assume without loss of generality that the fj have compact support and hence fj G 
L n (M. k ), for 1 < j < n. By (11071) with pj = n, we then have 

Wfl " • ' fn ~ h " " /n||l 

< ll(/r-/l)/|---/n!ll +H/l(/f-/2)/f •••» + •--+ HA' ■•/«-!(/: 

< ll(/l-A)l|n||/|||n---||^||n + ••• + || /l ||„ " ' ' II /n-1 1|» || Un ~ fn) \\n - 0. 

□ 

Part (i) of Proposition 13.141 is an immediate consequence of Lemma \3. 191 Next we show part 
(iii): By Lemma [3.191 and the independence of the s^, we have 

Kg A (5.KS)* -1 = E(Z s mMZ s m)aE(S £ (Z s mMZ s m))a---AE(S £ (Z s mMZ s m)) 
= E([Z s n B Z sf ] A El Z s *) A ■ ■ ■ A E Z sf )) . 
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Therefore, for a test form ip G V m ~ k,m ~ k (Q), we have 



k£a(s b k£)*-\ v 



/ (\Z s n B Z sf ] A S e (Z s * H Z sf ) A • ■ • A S £ (Z sf IE Z sf ), V 

JH°(M,L N ) k V 1 ' " " 



Hd 1N (s?) 



108) 



By Lemma [3.211 the integrand in (11081) is uniformly bounded, and hence by (1791) . Lemma [3.201 
and Lebesgue dominated convergence, we have 



(k£ A (^K^)* -1 ^) -> f fo* x XZ.x w 

X 7 JH°(M,L N ) k V 



n^( S f) 

b=i 



(K&,¥>) 



as e — >• 0, verifying part (iii). 

To complete the proof of the proposition, we recall from Proposition 13.101 that the current 
has L^ 1 coefficients (if m > 2) and hence by Lemma [5.221 



(S e K») j - (Kgy\n in L\Q), for 1 < j < m - 1 . 



^2i) \n 

Since L 1 convergence implies weak convergence, it follows from (iii) and (11091) that K^in 
(K^iYln and hence (ii) holds. This completes the proof of Proposition 13.141 



(109) 



-JVl 



□ 



3.5. Completion of the proof of Theorem 13.111 We recall that 

Vo\ 2m -2k[Z s N^ !S N nU] = (Z s n^ >s n, xu®k) = (Z s n^ jS n, Xu^k) as. (110) 
(To verify the second equality in (1 11 01) . we note that (Z s n n, Xdu&k) = almost surely, since 

1 '***' k 

E(Z s iv j s N,Xdu&k) = by Corollary 12.31 ) We now approximate Xu °y a sequence of C°° 
functions Xn '■ M — > R, n — 1, 2, 3, . . . , satisfying: 

• < X n < 1, 

• sup \dxn\ = 0(n), 

• Xn\u = 1, 

• Xn( w ) — for dist(f7, w) > 1/n. 

To construct Xn, we choose p G C°°(1R) such that p{t) — 1 for t < |, p{t) — for t > |, and 
< p < 1. Let Xnl^) = P( n dist(C/, w)). If <9£7 is smooth, then x n * s smooth, for n sufficiently 
large, and we can take Xn = Xn- Otherwise, the Lipschitz constant of x n is 0(n), and we can 
smooth x n to obtain our desired C°° function Xn- 

Then Xn — ► Xu pointwise, and hence for all (s± , . . . ,s^), we have by Lebesgue dominated 
convergence, 

{ z s?,...,s» iXn^fc) -> (^f,..., s f ,Xc/^fe) = Vol 2m - 2 fc[^iv,..., s A r n t/] as n -> oo . 
Therefore (again by dominated convergence), 

Var(Z sf i ___ iS jv, Xn<$>k) -> VarfVola™-^^,...,^ H Z7]) = Var(Vol 2m _ 2 fc[^f ,..., s f n f 7 ]) ( m ) 
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as n — > oo. To complete the proof of Theorem 13.111 it suffices by ( 11 111) and Theorem 13. 131 with 
if = Xn&k to show that 



B 1 B 2 Q N A {dAd^QNY A (EZ sN M EZ sN ) K 3 A (d[ Xn $ k } B 

<9l<9 2 Qiv A (^^^^g^)^ 1 A (EZ sN B EZ s ^- j A ($ fc IE $,) A (rf X n B 

:M 

-> - / d 1 d 2 Q N A {d 1 d 1 d 2 d 2 Q N ) j ~ 1 A (EZ s iv B EZ s iv) fc ^' A ($* B $ fc ). (112) 
Jauxdu 



'auxau 
To verify (TTT21) . let 

/ = BAQn A (didid&Qxy- 1 A (EZ s .v B EZ s iv) fc ^' A ($ fe B $ fc ) 
We must show that f\auxdu is an d 



/ A E rf X n) - / / ■ (H3) 

MxM JdUxdU 

By Lemma [3.71 we have 

|/(;z,™) | = O (iV J 'dist(z,w)- 2j+2 ) < O (iV J 'dist(z,w)~ 2m+2 ) . (114) 

Since <9£7 is a finite union of C 2 submanifolds of M of real dimension 2m — 1, it follows from 
(HUD that / is L 1 on <9t/ x dU . 

Let 5 > and consider the cut-off function Xs(z,w) = p(8~ 1 dist(z,w)), where p is as above. 
Then X s G C°°(M x M) for 5 sufficiently small, Xg(z, w) = if dist(z,w) > 5, and A^z, it?) = 1 
if dist(z, w) < 5/3. We decompose the integral in (1113ft : 

/ A (d X n H rf X n) = / A<S / A (d X n H d*n) + / (1 - A*) / A (t^n B rf X n) • (115) 

MxM J MxM J MxM 

Since (1 — Xg) f is smooth and Xn ^ Xu, it follows that 

/ (1 - \ S ) f A (rfXn B d* n ) = B d X n, (1 " X S ) f) 

J MxM 

-> (d Xf/ B dXt^, (1-X s )f) = - [ (1-X s )f. (116) 

Jauxdu 

(The minus sign in (jTTgjl is due to the fact that (ilB^B^) = (-l) dc g fi dc s^(A, 

and hence dxu E c?xc/ = [ — dU] B [—dU] = —[dU x dU], where [dU] denotes the current of 

integration over dU.) 

To complete the proof of (1113ft . we must show that the Xg f integrals are uniformly small. 

For Zq G M, n G Z + , 5 > 0, we write 

V(zo,n, 5) := {w G M : dist(^,w) < 5, w G Supp (rfx n )} 
C {u> : dist(z, w) < 5, dist(i7, w) < 1/ri}. 

Since <9?7 is piecewise smooth, we can choose 5q > 0, no G Z + such that for all zq G M: 

• the exponential map exp 2[) : T Zo (M) — > M is injective on the 5 -ball Bs (zo) — {v G 
•/•,(.!/) : |v| < 5 }; 
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• there exists real hyperplanes Pi, ... , P q , such that 



V(z , n, S ) C |J exp 2[) ({v + tuj G B So (z ) : v G Pj, \tj\ < 2/n}) 



;ii7) 



for all n > n , where Uj is a unit normal to Pj. 

Here, g is the maximal number of facets of the polyhedral cones locally diffeomorphic to open 
sets of dU, as described after the statement of Theorem ll.il (If dll is smooth, then q = 1.) 
Since j < m and \d\n\ = 0(l/n), we then have by (11141) and fl!17j) . 



\s(z , w)f(z , w) A dxn(w) 



M 



< 



f{z ,w) A dxn(w) 



{z }x V(z ,n,5) 



< Cn I dist(z ,w)- 2m+2 dVol 

'V(z ,n,S) 



M 



< Cn 



\x\ 



-2m+2 



dx 



{x<m. 2m :\x\<&, |xi|<2/n} 

2m + 2 dy = C"8 



< 4C" / \y 

'{y&R 2m - 1 :\y\<S} 



where C,C,C" are constants independent of z (but depending on m,U,N). Here, f(z ,w) is 
regarded as a (2m — l)-form (in the w variable) with values in T* 2m_1 (M). Therefore, 



\ S f A (d X n B d X n) 



MxM 



dXn(z) / X s (z,w) f(z,w) Adxn(w) 

{z£M:dist(U,z)<l/m} J{z}xM 



< C"5 / \d X n(z)\dVol du (z) 

J {z£M:dist(U,z)<l/m} 

< C"5 sup \dxn\ Vo\({z G M : dist(C7, *) < l/n}) . 

Since sup \dx n \ = 0(n) and the volume of the shell {z G M : dist(C7, z) < l/n} is 0(l/n), it 
follows that 



A 5 / A {dXn H 0?Xn) 



MxM 



< C'"5 V n > n . 



lis) 



Then (11131) follows from (j!15p . (I116p and (IllSp . which completes the proof of Theorem 13. Ill □ 



4. Variance of zeros in a domain: Proof of Theorems 11.11 and 11.41 

We now use Theorem 13.111 together with the asymptotics of the pluri-bipotential Qn to prove 
Theorem 11.41 

By Theorem 13.111 and Proposition 12.51 we have 



Var(Vol 2m _ 2fe [Z sf) ... )S N nU])=^2 (J) V f 



(U) , (119) 
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where 



if (CO 



2k-2j 



dUxdU 



-did 2 Q N {z, w) A [ - d 1 did 2 d 2 Q N {z, w)] J 1 



A 
1 



A' 



(m — /c)! 2 \ 7r 
where <3?fc is given by (|83l) . and 



A 



Tf A 



A $ fc (z) A $ k (w) 



du 



Tf(z) 



;.\ ( :. ' )| J 1 A 



(120) 



OJ (w 



-did 2 Q N (z,w) A [ - didid 2 d 2 Qi 

'{z}xdU 

e T* z j - l > j (M) . (121) 

(In (I12ip and below, we regard the integrand as an (m — 1, m)-form in the w variable with values 
in T* j -^(M) by identifying vr*^ A vr*^ e T*™ + ^ 2 ' m +^'(M x M) with y2 <g> ip G T*J- l ^(M) <g> 
T* m " 1 ' m (M), for G T*^(M), G T^ m ~ 1 - m (M) .) 
By Lemma 



B 1 B 2 Q N {z, w) A [«9 1 9 1 9 2 a 2 Q w (2, w)] j 1 = 0(A^ m ) , for dist(z, w) > b^jf , (122) 
where we choose b = i/2m + 3. Thus we can approximate Tj 1 (z) by restricting the integration 



in to the set of w G 9Z7 with dist(z,w) < V 1 ^- 

' 3 



To evaluate T^(z ) at a fixed point 2 G dU, we choose normal holomorphic coordinates 



{wi, . . . ,w m } centered at z and defined in a neighborhood V of z , and we make the change 
of variables Wj = ^= as in §2.21 Since uj = |<9<91oga = ~dd [\w\ 2 + 0(\w\ 3 )}, we note that 

1 i 
: dvj A dv k = -^^dd\v | 2 + O 



U[Zo + 



+ O 



N 



2N 



AT3/2 J 



(123) 



for |t>| < b^\og N. We then have 
Tf(zo) = W- 



-did 2 Q N (z Q ,z Q + J=) 



{\v\<bVWN:z +^=£du} 

A [ - d 1 B 1 d 2 d 2 Q N (z , z + ^Y" 1 A [(!&9M 2 )^ + O (£)] 
Applying the asymptotics of Lemmas I3.8H3.9I to (I124p . we obtain the formula 



(124) 



A (Var^y- 1 A ( % -dd\v\ 2 ) + 0(Ar x / 2+£ ) 



m-j 



(125) 



We first consider the case where dU is C 2 smooth (without corners). We can choose our 
holomorphic normal coordinates {wj} so that the real hyperplane {Imwi = 0} is tangent to 
dU at z . We can then write (after shrinking the neighborhood V if necessary), 

U f)V = {w G V : Im W\ + <p(w) > 0} , 
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where ip : V — > R is a C 2 function of (Re w±, w 2 , ■ ■ ■ , w m ) such that <^(0) = 0, dip(0) = 0. 
We consider the nonholomorphic variables 

w = t(w) := (wx + iip(w),w 2 , . . . , w m ) , (126) 

so that dU = {Imwi = 0}. We next make the change of variables 

Nw =v 1+0 



v = t n (v) := V N r 
in the integral (I125p to obtain 

Tf(zo) = N 2j ~ m - x ' 2 



lF"{±\v\ 2 )8{z-v) Ad\v\ 2 



A (VaC(tO) i-1 A (idd\v\ 2 ) m ~ J + 0(AT-i/2+ £ ) 



(127) 



where 



G R x C" 1 " 1 : \v\ <{b- l)v/logiv} C C |t> G R X C^ 1 : |v| < {b + l)y/]ogN^ 

By (JTOD and (JTSD, we have F''(±H 2 )H 2 {Var£(v)p'- 1 = 0(e~H 2 ) for |v| > 1, and hence 



[«l>(6-i)VEgW 



F"(||£| 2 ) 9(2 • tT) A <9|£| 2 A (Var£(tO) i-1 A (i<9<9|v| 2 ) m ~ J 

= O (jV^M = O (N- 1 ) . 

Thus we can replace the B 2 ^ 1-1 integral in (1127ft with the affine integral over R x C m_1 , so that 

Tf (z ) = iV^— 1 ^ [Tf (^o) + 0{N- l ' 2+e )\ , (128) 

where 



Tr(^o) 



1 F"(|M 2 ) <9(z • v) A <9M 2 A (Var^V 1 A (^dd\v\ 2 ) 



\2\ m -3 



4 V2 



- d(z ■ v) A B\v\ 2 A (Var^y'" 1 A (f<9<9M 2 ) m 3 . (129) 



47r 2 / , pW 



Since only the last two terms of Var^ contain neither d(z-v) nor d\v\ 2 , formula ( 1129f) simplifies 
to: 



T3°(*b) 



Thus, 



-1) 



(4vr 



2V 



1 



, — A did 2 

A [dd{z ■ v) A dd{z ■ v) + 99|^| 2 A ddlvl 2 ] 1 ' 1 A (f<9<9M 2 ) m ~ J . (130) 



(T°° A u rn j )(z ) = c mj dxi A (§tte 2 A dz 2 ) A ■ • ■ A (§Gfe m A dz m ) = c mj dVo\ aUtZo , (131) 

where c m j is a universal constant, which we compute in §4. II below. 
Substituting (11281) and (TTrTTT) in (jT20l) . we have 



if (to 



1 



(m — /c)! 2 \7r 2fc 2j 



yy2fc-TO-l/2 



/ c mi t/VoW. 20 + O (iV- 1 ^) . (132) 
vac/ 



NUMBER VARIANCE OF RANDOM ZEROS ON COMPLEX MANIFOLDS 



37 



Combining ( I119p and ( 11321) . we obtain the formula of Theorem 11.41 with 

k 



(m — k)\ 



J J 7T 2k 2 ^ 



(133) 



for the case where dU is smooth. 

We now verify the general case where dU is piecewise smooth (without cusps) . Let S denote 
the set of singular points ('corners') of dU, and let Sn be the small neighborhood of S given 
by 



'N 



zedU : dist(z,£) < 



b'^\ogN 



N 



where b' > is to be chosen below. We shall show that: 

i) (I128p holds uniformly for z Q G dU \ Sn] 

ii) sup \Tf(z)\ = O (jv 2 -?-™- 1 ^) , for 1 < j < k. 



zedu 



Let us assume (i)-(ii) for now. Since Vol 2m _iSjv = O ( ~^p^ j , the estimate (ii) implies that 

0(N 2j - m - l+e ) , 



and hence by (11201) . 



V N (U) 



(m — k)\ 2 \ 7T 



fN \ 2k ~ 2j 



dU^S N 



Tf A 



10 



O (N 



2fc-m-l+e 



)■ 



It then follows from (i) and (11311) that 

yy2fc— m-l/2 



V»{U) 



2 ~2k-2j 



T°° A u m ~ j + 0{N~ l/2+£ ) 



au^s N 



(m — k)\ 2 7r 2k 2 i 
Then by ( TTT91 

Var(Vol 2m _ 2fc [Z sfi ... jSf nf/]) = N 2k - m ~ 1 ' 2 

7y2fc-m-l/2 



(m — k)\ 2 7T 
_ T\j2k-m-l/2 

[Vol(9[/x^) + 0(iV- 1 / 2+£ )] . 



^VoL.^f/ \ S N ) + 0(N—2 
v mk Vo\ 2m ^{dU)+0{N-^ 2£ 



which is our desired formula. 

It remains to prove (i)-(ii). To verify (i), for each point zq G dU \ S, we choose holomor- 
phic coordinates {wj} and non-holomorphic coordinates {wj} as above. We can choose these 
coordinates on a geodesic ball V ZQ about z of a fixed radius R > independent of the point z Q , 
but if zq is near a corner, dU will coincide with {Imroi = 0} only in a small neighborhood of 
z . To be precise, we let D zo denote the connected component of V zo D dU \ S* containing z . 
Then we choose G C 2 (V Z0 ) with 9?(0) = 0, dip(0) = 0, such that 

{w G Ko : Imwi + (p(w) = 0} = {Imij = 0} D -D 2Q . (134) 

We let 

dist(2,5) 

C = SUp ; r > 1 . 

zeau^s dist(z,dE/ \ D z ) 
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Choose JV > such that bJ ^jM < r ; then 

G dU : dist(z , w) < C V zo , for N > N . 

We recall that our assumption that dU is piecewise C 2 without cusps means that U is locally 
C 2 diffeomorphic to a polyhedral cone, which implies that C < +oo. We now let b' = Cb, where 
b = \j2m + 3 as before. 

Consider any point z G dU \ Sn, N > N . Then 

a-+( Pitt n w dist(^ fo^ogW feylojiV 

Thus by our far-off- diagonal decay estimate (11221) . the points in dU \ D ZQ contribute negligibly 
to the integral in (11251) . so that integral can be taken over the set 



jM < by/log N : z + -^=eD Z0 \ 



which is mapped by into M. x C m 1 . Then (11271) holds, and (11281) follows as before. 
To verify (ii), we must show that the integral in the right side of (I125p . 



Tf(zo) := / \ F"(\\v\ 2 ) B(z ■ v) A B\v\ 2 A (Var^" 1 A (f dd\v\ 2 ) m J 

J{\v\<bVWN:z + -^=(idu) 

is 0(N £ ) uniformly for zq G dU . By Lemma [3.9[ Var^(f) = 0(\v |~ 2 ). Furthermore, 



1 



(v ■ dz) A (v ■ dv) 



m \v 



2 



;M a - 1 

(using Euclidean norms in the z and v variables ), and hence 



< — j — f« < \pm 

- P \v\ 2 _ 1 - v 



|Tf (z )| < A jm \v\-V+ 2 dVol^v) , (135) 

for universal constants A, m , where Vol B denotes Euclidean volume. Rewriting (I135P in terms of 
the original variables w = z + we have 



|Tf (z )\ < A jm N-^l 2 f \w- z \-^+ 2 dVol*^ 



For each point P G <9?7, we choose a closed neighborhood Vp G M of P and a C 2 diffeomor- 
phism pp : Vp —> IR 2m mapping Vpfldf/ to the boundary of a polyhedral cone Kp C IR 2m . Then 

for iV sufficiently large, for all ;z G <9^7 5 the set |u» G <9£7 : \w — z \ < byj^j^- \ is contained in 

one of the Vp. We then make the (nonholomorphic) coordinate change w = pp(w). Since the 
diffeomorphisms pp have bounded distortion, we then have 



|Tf (z )\ < A> m AT-^-+V2 / |~ _ ~ o| -2 i+ 2 JVorJ^) . (136) 

J\wedK P -.\w-zo\<b^^Lj 
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Let q G Z + be the maximum number of facets in OK p. We easily see that 
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J < w 



{wedK P -.\w-z \<b^f} 



\W ~ Zq\ 



-V +2 <No^ m _ x {w) < q [ \x\~ 2 ^ +2 dx 

const, x I -j— J . (137) 



Combining (j!36p — (j!37j) . we conclude that Tj^ (20) = 0(N £ ) and thus (ii) holds, which verifies 
the formula of Theorem 11.41 for the general case where dll has corners. 



4.1. Positivity of the constants v m k- To complete the proof of Theorems 11.11 and 11.41 we 

must show that the leading coefficients v m ^ are positive. In fact, we shall show that each of the 
coefficients c m j (defined by (11311) ) is positive, and then the positivity of the v m ^ follows from 



We begin by computing the coefficient u m \ in the co dimension-one case: By (11291) . 



T ~(-> = ^E 

(m - 1)! 



xC" 



m—l 



d Zj 



4tt 2 



E 

3=1 



VjVi 



xC r - 



-1 eM 2 



dz 



3 ■ 



(138) 



By (!TT9j) - (pD]l with k — l and ffT28l) . we have 

N 3/2-rn 



Var(Vol 2m _ 2fc [Z sf nU]) 



(m-1)! 2 [J 9U 



Tf (z) A ujiyz)™- 1 + O {N- l ' 2+£ ) 



(139) 



Since cfVobjc/^o) = dxi A ( . m j_ 1 ^, u;(zo) ?n only the j = 1 term in (I138p contributes to the integral 
in (11391) . and we then have 



Var(Vol 2w _ 2 *[Z sf fW]) = iV 3 / 2 "™ / v ml dVol 9u + O (N-^ 2+£ ) 

Uau 



(140) 



where 



1 



1 



4tt 2 ./ R 2m-i eM 2 - 1 
1 



2 7J -m-l/2 roc ^2m 



4vr 2 (2m - 1) T(m - 1 /2) J e r2 - 1 

m— 5/2 00 /-00 



47T 2 (2m - 1) J R2m -i eM 2 - 1 
dr 



dv 



7T 



4T 



(m+l/2)^7 



7T 



m-5/2 



E 



T(m + l/2) 7T 



m-5/2 



4T(m+l/2)^ 2A; m + 1 /2 
as stated in the theorem. 



C(m+| 
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We now determine the c m j, for 1 < j < m. By ni3Q[) ( TT3Tj) . 

c m j dxi A (^dz 2 A dz 2 ) A ■ ■ ■ A [\dz m A dz m ) = c m j dVolgu 
(-l)i- 1 



A [{dz ■ dv) A [dz ■ dv) + (dz • dz) A ((dv • dv)Y~ l A [§00|*| 2 A §<9<9M 2 ] m ~\ (141) 

where t> ■ dz = v a dz a , dz ■ dv = ^ dz a A dv a , etc. As in the codimension one case, only the 
term v\dz\ A dv\ in (v ■ dz) A (v ■ dv) contributes to the integral, and we obtain 

(. \ 2m— 2 ~[ r 2 

\) L„_, w^v d ' z ' A *■ A (A + Brl A B ""' ' (142) 



where 



01 



A = , Aq,^ = dz a A d,zg A dvp A (if 



B = _B a/ g , B a p = dz a A dz Q A dv/3 A du/3 . 

a,/3=2 



Writing 



j'-i / • _ i\ 

(A + S) 3 ' -1 A 5 m ^' = ^ ( J ) ^ A B™-'- 1 , (143) 
we see that it suffices to show that (f ) 2 ™ 2 A £ A B" 1 ^^ 1 is positive for < £ < m — 1. We have 



4* AjB m-/-i = ^ ... ^ ^A-Ai^AB^A-AB^^. (144) 

ai,/3l=2 a m _i,/3 m _i=2 

We claim that each term in (11441) equals B22 A B 33 A ■ • • A B mm if ai, . . . , a m _i is a permutation 
of 2, . . . , m, and . . . , /3 m _i is obtained from a±, . . . , a m _i by permuting ol\, . . . , and also 
permuting ae+i, . . . , a m -i- Otherwise, the term clearly vanishes. To verify the claim, by simul- 
taneously permuting the coordinates {z a } and {v a }, we can assume without loss of generality 
that a p = p + 1 (1 < p < m — 1); starting with (3 P = a p , where the claim is a tautology, we 
note that a transposition in f3\, . . . , transposes two (dvp A dz^s, while a transposition in 
/3i + i, . . . , Pm-i transposes two (dvp A dvp)\ leaving the sign unchanged. 
If follows that 

(|) 2m_2 dzt A dvt A A e A B m ' 1 - 1 = (m - (m - 1 - 1)! dVoW A dVol RxC m-i, 

XjL/ ai/x(RxC m - 1 ) " 

(145) 

and thus the c mj - are positive, completing the proof of Theorem 11.41 □ 
Combining f)133j) . (j!42p . H143j) and (I145p . we obtain the explicit formula: 



n m-2k-i/2 k } ( m _ -q j * n(m,j) Z" 00 r 2m dr 



4T 



(m + 1/2) (m - A;)! 2 ^ j (A; - j)\ J (e r2 - 1) 



2-^ nth - n\\ L _ ' 
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where n(m,j) G Z + is given by 



(m — 



= E^r- < 147 > 



5. Appendix: Proof of Theorem 12.41 

In this appendix, we sketch the proof of the off-diagonal Szego asymptotics theorem. The 
argument is essentially contained in [SZ2] , but we add some details relevant to the estimates in 
Theorem 12.41 

The Szego kernels Un(x, y) are the Fourier coefficients of the total Szego projector U(x, y) : 
£ 2 {X) -> H 2 {X); i.e. U N (x,y) = ± J e~ iN0 U{e i9 x,y) dB. The estimates for U N {z,w) are then 
based on the Boutet de Monvel-Sjostrand construction of an oscillatory integral parametrix for 
the Szego kernel: 

II(z,y) = S(x,y) + E(x,y) , 

(148) 

with S{x,y)= f °° e it1p( - x ^s{x,y,t)dt, E(x, y) G C°°(X x X) . 

The amplitude has the form s ~ XlfcLo t m ~ k Sk{x, y) G S m (X x X x The phase function ip 
is of positive type, and as described in |BSZ2j . is given by: 



if)(z, 0, w,ip) = i 



^a(z)s/a(w) 



(149) 



where a G C°°(M x M) is an almost holomorphic extension of the function a(z,z) := a(z) on 
the anti-diagonal A = {(z,z) : z G M}, i.e., da vanishes to infinite order along A. We recall 
from f[T5]) that a(z) describes the Hermitian metric on L in our preferred holomorphic frame at 
z , so by fl2~4"l) . we have a(u) = 1 + \u\ 2 + 0(\u\ 3 ), and hence 

a(u,v) = 1 + u ■ v + 0(\u\ 3 + \v\ 3 ) . (150) 

For further background and notation on complex Fourier integral operators we refer to [BSZ2J 
and to the original paper of Boutet de Monvel and Sjostrand [BS] . 

As above, denote the iV-th Fourier coefficient of these operators relative to the S* 1 action by 
II/v = Sn + En- Since E is smooth, we have E^(x,y) = 0(N~°°), where 0(N~°°) denotes 
a quantity which is uniformly 0(N~ k ) on X x X for all positive k. Then, En(z,w) trivially 
satisfies the remainder estimates in Theorem 12. 4[ 

Hence it is only necessary to verify that the oscillatory integral 

f'I'K />oo t>2n 

S N (x,y)= / e- im S(e i0 x,y)d6= / / e- im+it ^ ex > y) s(e ie x,y,t)d6dt (151) 



Jo Jo Jo 

satisfies Theorem 12.41 This follows from an analysis of the stationary phase method and re- 
mainder estimate for the rescaled parametrix 



2tt 



^ -=,0;-=,0 =JV / / e V ^ ^>) s [-=,9; -7=,0,JVt ) dOdt , 

N \Vn Wn J Jo J \Vn Wn ) 

(152) 

where we changed variables t i— > Nt. For background on the stationary phase method when 
the phase is complex we refer to [Ho] . We are particularly interested in the dependence of 



42 



BERNARD SHIFFMAN AND STEVE ZELDITCH 



the stationary phase expansion and remainder estimate on the parameters (u, v) satisfying the 
constraints in (i)-(ii) of Theorem 12.41 

To clarify the constraints, we recall from [SZ2, (95)] that the Szego kernel satisfies the fol- 
lowing far from diagonal estimates: 

\V 3 h Ii N {z,w) \=0{N- K ) for all j, K when dist > 



N 



Hence we may assume from now on that z = zq + w = zq + ^7= with 

|u| + \v\ < 5N 1/6 

for a sufficiently small constant 5 > 0. 

By (1149D — (I150p . the rescaled phase in (I152p has the form: 



(153) 



(154) 



* := tip 



11 e-^n 



N 



N 



9 = it 



JO 



(155) 



it 

—ip 2 (u,v)e v 



tffl 



u 



V 



N VN 



(156) 



and the TV-expansion 

* = - e i9 } 
where 

tp 2 {u,v) — u • v — \{\u\ 2 + |f | 2 ) = — ||w — v\ 2 + z Im {u ■ v) . 
After multiplying by iN, we move the last two terms of (11561) into the amplitude. Indeed, we 
absorb all of exp{(^2 + iNR^)te t9 } into the amplitude so that (11521) is an oscillatory integral 

poo r2ir 

N 



JO 



with phase 



and with amplitude 



V(t,6) := it{l - e ld ) - 9 



A(t,9; Zq,u,v) := e * K ' ' 3^'v^v ; s( 



u e V 



0,Nt). 



(157) 



158) 



(159) 



N VN' 

The phase ^ is independent of the parameters (u, v), satisfies Re (i^>) = —t(l — cos#) < 
and has a unique critical point at {t — 1, 9 — 0} where it vanishes. 

The factor e ie ' ^ 2 ( u < v *> is of exponential growth in some regions. However, since it is a rescaling 
of a complex phase of positive type, the complex phase iN^f + te ip 2 (u, v) is of positive type, 

Re (iNV + te ie ip 2 (u, v)) < (160) 

once the cubic remainder Nte l6 R^(^, is smaller than iN^/ + te l9 ip 2 (u,v), which occurs 
for all (t, 6, u, v) when (u, v) satisfy (11541) with 5 sufficiently small. 

To estimate the joint rate of decay in (N, u, v), we follow the stationary phase expansion and 
remainder estimate in Theorem 7.7.5 of |Hoj . with extra attention to the unbounded parameter 
u. 

The first step is to use a smooth partition of unity {pi(t, 9), p 2 (t, 9)} to decompose the 
integral (11521) into a region (1 — e, 1 + s) t x (— e, e)$ containing the critical point and one over 
the complementary set containing no critical point. We claim that the p 2 integral is of order 
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N~°° and can be neglected. This follows by repeated partial integration as in the standard 
proof together with the fact that the exponential factors in (11601) decay, so that the estimates 
are integrable and uniform in u. 

We then apply [H6J Theorem 7.7.5 to the pi integral. The first term of the stationary phase 
expansion equals ]S[ m e teXe ^{u,v) an( j remainder satisfies 

\Rj(P ,u,v,N)\ <CN~ m+J suv\D? )ePl A(t,9;P ,u,v)\. (161) 

H<2J+2 tfi 

From the formula in (11591) and the fact that s is a symbol, A has a polyhomogeneous expansion 
of the form 

K 



A(t,0;P o ,u,v) = Pl {t,6)e teWMu ' v) N m 



N- n/2 f n (u } v; t, 6, P ) + R K (u } v, t, 6) 

.71=0 

K+l , 



\V j R Nk (u,v)\ < C ]k£b e^ u \ +H ^N—t-). (162) 

The exponential remainder factor e e( -' ,1 ' 2+ ' , '' 2 - ) comes from the fact Ree* 6 -?/^ = cos#Re^ — 
smOlmip with Keip < and |sin#| < £ on the support of p\. Hence, the supremum of 
the amplitude in a neighborhood of the stationary phase set (in the support of pi) is bounded 
by e £ l Im ^ 2 L The remainder term is smaller than the main term asymptotically as iV — > oo as 
long as (u, v) satisfies (I154p . Part(i) of Theorem 12.41 is an immediate consequence of (I162p since 

e e(\u\ 2 + \v\ 2 ) < N e for | w | + | v | < V 4 ^7\T. 

To prove part (ii), we may assume from (I153p - (I154I) that \/\og N < \u\ + \v\ < 5 N 1 ^ . In this 
range the asymptotics fl 1 6 2 [) are valid. We first rewrite the horizontal z-derivatives Jj- as Uj 

derivatives, which for L N have the form \fN t£- — NAj(^=) and thus contributes a factor 

of y/N \ We thus obtain an asymptotic expansion and remainder for V^IIjv^, w) by applying 
to the expansion (i) with k = 0: 

nf («, 6- v, cp) [1 + N-^Rmtu, v)] . 

The operator contributes a factor of N^ 2 to each term, and thus 

\V j h n N (z,w)\ = O (N m+j / 2 e'^^ 

= 0(N~ k ) uniformly for \u\ 2 + \v\ 2 > (j + 2k + 2m + e')\og N , 
where s' = (j + 2k + 2m+ l)e. □ 
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